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HYMO? 

PROBLEM-ORIENTED COMPUTER LANGUAGE 
FOR HYDROLOGIC MODELING 

Users Manual 

By Jimmy R. Williams^ and Roy W. Hann» Jr.^ 


INTRODUCTION 


HYMO ( 7 )^ is a problem-oriented computer 
language for modeling surface runoff and sedi- 
ment yield from watersheds. The language is 
called HYMO from the words “hydrologic 
model.” HYMO was designed for planning flood 
prevention projects, forecasting floods, and re- 
search studies. It consists of a main program and 
16 subroutines written in FORTRAN IV, but it 
can be used by hydrologists with little knowl- 
edge of computer programing. The language 
provides 17 commands for the hydrologist to 
use in any sequence for application to any 
watershed. 

HYMO was designed to transform rainfall 
data into runoff hydrographs and to route these 
hydrographs through streams and valleys or 
reservoirs. It will also compute the amount of 
sediment produced by a storm at any point on a 
watershed. It will be useful to research hydrolo- 
gists in studying the effects of watershed and 
storm characteristics on the flood hydrograph. 
HYMO is also a good research tool for testing 


hydrologic procedures; for example, a new 
flood-routing method could be added to HYMO 
and tested easily, because the inflow hydro- 
graphs and the rating curves are available in a 
HYMO program. 

HYMO is flexible. Present hydrologic proce- 
dures can be modified or deleted, and other 
hydrologic procedures can be added by hydrolo- 
gists familiar with FORTRAN IV programing. 
Adding a new command simply requires the 
addition of a new subroutine. 

HYMO is efficient, practical, and generally 
applicable. HYMO programs can be written and 
the results interpreted by hydrologists who have 
no conventional computer programing experi- 
ence. The hydrologic procedures used in HYMO 
are practical — required inputs are easy to obtain 
for most watersheds. 

HYMO was written for the IBM 360-65 
computer, but it could be run on an IBM 1130 
with little modification. The storage require- 
ment is about 73 K. 


OPERATION OF HYMO 


HYMO consists of a main program and 16 
subroutines. The HYMO card deck is set up in 
the following order: 


'Hydraulic engineer, Coilege Station Area, Southern Region, 
Agricultural Research Service, U. S. Department of Agri- 
culture, Riesel, Tex. 

"Head, Environmental Engineering Division, Civil Engineer- 
ing, Texas A&M University, College Stotion, Tex. 

"Italicized number in parentheses refer to items in "Liter- 
oture Cited" preceding the appendix. 


1. Main progrsim. 

2. Subroutines. 

3. A data card containing the number of 
commands in the command table. 

4. A data, card containing the ZALFA array. 

B. Seventeen data cards containing the com- 
mand table. 

6. The users program deck consisting of 
program and data cards. 

A printout of the main program, subroutines, 
ZALFA array, and command table is given in 
the appendix, 
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The main program reads the command table gram data card with the Z ALF A array . The data 
and then calls the HONDO subroutine (2) to are placed in an array and returned to the main 
read a program data card. Subroutine HONDO program. Based on the command number, the 
determines the command name and number by main program calls the proper subroutine to do 
comparing columns 1 through 20 of the program the desired calculations. When the calculations 
data card with the command table. Then are complete, control is returned to the main 
HONDO determines individual data items by program, and HONDO is called again to read the 
comparing columns 21 through 80 of the pro- next program card. 


HYDROLOGIC PROCEDURES USED IN HYMO 


The procedures used in HYMO were selected 
because of their accuracy, general applicability, 
practicality of inputs, and computational effi- 
ciency. For most watersheds the input is easy to 
obtain, and the procedures produce reasonably 
accurate results without excessive computer 
time. 

Hydrograph Computation 

When flood routing is performed, a watershed 
is divided into many small areas according to its 
hydraulic characteristics. The hydrographs from 
these areas must be estimated, since streamflow 
measurements are seldom available. A procedure 
for computing unit hydrographs was developed 
previously (4). A modification of this procedure 
is used in HYMO. Unit hydrographs are divided 
into three parts for computation (fig. 1). From 
the beginning of rise to the inflection point, /q > 
the hydrograph is computed by the two- 
parameter gamma distribution equation 



Figure 1. Dimensionless unit hydrograph. 


where q = flow rate in cubic feet per second at time (, 
Qp “ peak flow rate in cubic feet per second, 
ip = time to peak in hours, 
and n = dimensionless parameter. 

From fo to t\ (fi = (q + 2K) the hydrograph 
is computed by the recession depletion equation 

K 

Q “ Qoe ( 2 ) 

where <7o = flow rate at the inflection point, 

Iq » time at the inflection point, 
and K = recession constant in hours, 

From fi to oD the recession depletion equa- 
tion becomes 

f|-f 

Ki 

^ Qie ( 3 ) 

where </| = flow rate at , 

and /T] » 3K = second recession constant. 

The dimensionless shape parameter, n, is a 
function of K/tp, as shown in figure 2. The peak 
flow rate is computed by the equation 

BAQ 

" tp (4) 

where Z? = a watershed parameter, a function of n as 
shown in figure 3, 

A - watershed area in square miles, 
and Q - volume of runoff in inches, 

Therefore, the entire unit hydrograph can be 
computed if K and tp are known. K and tp 
can be determined by hydrograph analysis (4) 
for gaged watersheds. To compute K and tp for 
ungaged watersheds, HYMO uses the equations 

K = 27.0^0.23l5ip-0.777('Z,/lv;0.124 /gj 
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and tp = 4.63/lO-4225Li^-0.46^Z,/H/;0.133 (6) 

where SLP = difference in elevation in feet, divided by 
flood-plain distance in miles, between 
watershed outlet and most distant point 
on the watershed, 

and LjW = watershed length-width ratio. 

Storm hydrographs are computed by con- 
volving unit hydrographs with incremental 
source runoff. To compute incremental source 
runoff, the mass rainfall curve is broken into 
equal time increments, and the Soil Conserva- 
tion Service (SCS) rainfall-runoff relationship 
(3) is applied. The SCS rainfall-runoff relation- 
ship is expressed in a set of numbered curves. 
The SCS National Engineering Handbook (J) 
provides detailed instructions for selecting the 
proper curve number. 

Hydrographs computed by this procedure 
compared closely with measured hydrographs 
from 34 watersheds located in Texas, Oklahoma, 
Arkansas, Louisiana, Mississippi, and Tennessee. 
The watershed areas ranged from 0.5 to 25 
square miles. 

Flood Routing 
Streams and valleys 

The variable storage coefficient (VSC) flood- 
routing method (5) was selected for HYMO. 



Figure 2. Relationship between dimensionless shape 
parameter and recession constant/time to peak. 



Figure 3. Relationship between dimensionless shape 
parameter n and watershed parameter B. 


The VSC method has been revised ((5) to 
account for the variation in water surface slope 
during a flood. The revised VSC method is about 
as accurate as the implicit method (/) and has 
the general applicability of simpler storage 
methods. Although an iterative solution is used, 
the VSC method requires little computer time 
and is free of convergence problems. 

The VSC routing equations are 


C2 - 
Cl = 


c,\!, ♦(i-i) ol 

2Af 

(7) 

(8) 

2 T 2 +^t 

2^t 

(9) 

2Ti+^t ■ 

f ^ ^ 


^^ISOOfFij + Vq^JJ 


I lkslpc, 

(10) 

-^0, / 

( ^ ] 



1800 (Vi^ + Vo^)) 
L v- SLPq 




SLPo +£>1, -X>02 


% 


( 11 ) 
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In these equations subscripts 1 and 2 refer to 
the beginning and end of the time interval a/; the 
units are cubic feet per second for flow, hours 
for time, feet per second for velocity, and feet 
for length and depth. The symbols are defined as 
follows: 


/ = inflow rate. 
0 = outflow rate. 


h = 


■y = average inflow rate. 


C = storage coefficient. 

T = travel time through the reach. 
L = reach length. 

V = velocity. 

SLPq - normal slope. 

D = depth. 


Since Ti and C 2 are dependent upon O 2 , an 
iterative technique is required to solve the 
routing equations. In equation 7, Iq and 0i are 
known, and Cj can be computed from equation 
9. This leaves only O 2 and C 2 as unknowns. 0\ 
can be used as a first approximation of O 2 . The 
normal depth and velocity for the approximate 
value of O 2 are entered into equation 11 for 
computing T^. Then equation 8 is used to 
compute C 2 . The second approximation of O 2 is 
then obtained from equation 7. This iterative 
process continues until the difference between 
successive O 2 values is acceptable, HYMO is set 
to accept differences of 0.1 percent or less. 
Usually about four iterations are required. 
Reservoirs 

HYMO uses the storage-indication method (5) 
to route floods through reservoirs. This method 
has been widely used and accepted because it is 
practical and accurate. The SCS National Engi- 
neering Handbook gives detailed instructions for 
using the method. 


Rating Curves 

Rating curves must be available at enough 
locations along a valley to adequately describe 
the hydraulics of the stream and valley. Most of 
these rating curves must be computed because 
there are never enough measured rating curves. 

HYMO uses Manning’s equation to compute 
the normal flow-rating curves that are used in 
the VSC flood-routing method. The normal 
flood-plain slope is determined for each valley 
section by plotting a profile of the flood plain. 
The normal channel slope is determined by 
plotting a profile of the flood plain with channel 
distances. 


Sediment Yields 

The universal soil loss equation (8) was 
modified to compute the sediment yield for 
individual storms on watersheds. The modified 
equation is 

S=9Q(qpXQ)^-^^ xKxCxPxLS, ( 12 ) 


where S 

% 

K 

C 

P 

and LS 


sediment yield in tons, 

peak flow rate in cubic feet per second, 

volume of runoff in acre-feet, 

the soil-erodibility factor, 

the cropping-management factor, 

the erosion control practice factor, 

the slope length and gradient factor. 


Detailed instructions for determining K, C, P, 
and LS are given by Wischmeier and Smith (5). 

Since equation 12 was designed to compute 
sediment yield from watersheds, a delivery ratio 
is not needed. The delivery ratio is built into 
equation 12 by including the peak flow rate. 
Many of the watershed characteristics that in- 
fluence the peak flow rate also affect the 
delivery ratio. Equation 12 has performed well 
under limited testing, but future refinements are 
expected. 


RULES FOR USING HYMO 


The reader should refer to the table, ‘‘Ex- 
ample Input for HYMO Commands,” as he 
follows the narrative description of the rules. 
The example HYMO program that is presented 
near the end of the manual may also be helpful. 

General Rules 

HYMO commands are expressed in the first 
20 columns of the data card, and columns 21 


• through 79 are used for numeric data and 
‘keywords. Column 80 is reserved for a page^ 
change code (an asterisk in column 80 causes the 
card to be printed on a new page). Continuation 
cards are allowed when 59 characters are insuf- 
ficient to express the data. 

The data can be written in any format, but at 
least one blank space must be left between data 
items. A decimal is required for numbers con- 
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taining fractions, but not for whole numbers, may be used at any point in a HYMO program 
Keywords can be written with the data to by punching an asterisk in column 1 and the 
describe individual data items. Comment cards comment in columns 2 through 79. 


Example input for HYMO commands 

Command Required input 

START RAINFALL BEGINS AT 12.5 HRS PUNCH CODE=l 

STORE HYD ID=1 HYD NO=301 DT=.2 HR DA=1.5SQMI 

FLOW RATES= 0 10 50 100 500 1000 1800 2000 1900 1500 1200 
1000 800 600 600 400 300 200 100 60 10 1 

COMPUTE HYD ID=2 HYD NO=302 DT=.6 HR DA=2.1SQMI CN=90 

HT=100PT L=3.3MI 

MASS RAINFALL = 0 .31 .61 1.04 1.84 2.74 3.06 3.46 4.33 4.76 
PRINT HYD ID=2 CODE=l 

PUNCH HYD ID=2 

PLOT HYD ID 1=3 ID 11=4 

ADD HYD ID=4 HYD NO»101 ID 1=6 ID 11=6 


STORE RATING CURVE ID-2 VS NO=16 


ELEV 

AREA 

FLOW 

496.6 

0 

0 

497 

2 

1 

498 

9 

19 

499 

19 

52 

600 

30 

98 


COMPUTE RATING CURVE ID=1 VS NO=10 NO SEGS=3 MIN ELEV=482 FT MAX ELE V=492 FT 

CH SLP=.006 FP SLP=.0076 N=.06 DIST=176 FT 
N=-.03- D1ST=20B FT N=.0B DIST=460 FT 


DIST 

ELEV 

DIST 

ELEV 

DIST 

ELEV 

DIST 

ELEV 

0 

492.0 

100 

490.0 

176 

484.0 

188 

482.0 

190 

482.0 

20B 

484,0 

250 

486.0 

276 

488,0 

310 

490.0 

4B0 

492.0 






COMPUTE TRAVEL TIME ID=3 REACH NO=8 NO VS=6 L=4600 FT SLP=.007B 
ROUTE ID=3 HYD NO=8 INFLOW ID=6 DT=.26 HR 


ROUTE RESERVOIR ID=6 HYD NO-BOl INFLOW ID=1 


OUTFLOW (CPS) 

STORAGE (AC FT) 

0 

0 

22 

633 

200 

656 

1000 

• 601 

2000 

648 

3000 

694 


ERROR ANALYSIS ID 1=3 ID 11=6 

SEDIMENT YIELD ID=B SOIL=.34 CROP-.B EP=.6 LS=.3 




Six hydrographs can be stored in a HYMO punched for use in the second routing. More 
program at a time. The hydrographs are identi- than one routing is usually required, 
tied by storage location numbers 1 through 6. Two commands, RECALL HYD and STORE 
Therefore, the same storage location number TRAVEL TIME, were designed to be computer 
must be used for many hydrographs in a HYMO punched for second routings; consequently, 
program. This is especially true when routing is these commands do not appear in the table, 
done through lai’ge watersheds. However, no The STORE HYD command is used to store 
more than six hydrographs are ever needed at the coordinates of a hydrograph in the corn- 
one time because HYMO programs begin at the puter. It can be used for storing measured 
head of a watershed and work downstream hydrographs or hydrographs computed by 


through one reach at a time. When a storage 
location number is used to store or compute 
another hydrograph, the first hydrograph is lost. 
The user should be sure that the hydrograph will 
not be referred to again before using the storage 
location number for another command. 

To store, compute, or route a hydrograph, the 
user must specify the time increment. There are 
no rigid rules about selecting the time incre- 
ment, but generally it should not be greater than 
one-fifth of the time to the peak of the 
hydrograph. This rule usually provides enough 
points to adequately define the hydrograph. All 
hydrographs are limited to 300 points. 

For the commands “STORE HYD,” “COM- 
PUTE HYD,” “ADD HYD,” “ROUTE,” and 
“ROUTE RESERVOIR,” the user must specify 
the number of the outflow hydrograph. The 
hydrograph identification numbers are used to 
designate specific routing reaches, incremental 
areas, reservoirs, and partial hydrographs. The 
partial hydrograph number is given to all hydro- 
graphs other than outflow hydrographs from 
reaches, incremental areas, or reservoirs. The 
identification numbers for each group are 


Reaches 1-100 

Partial hydrographs 101-300 

Incremental areas 301-500 

Reservoirs 501+ 


Command Rules 

The first command for any watershed is 
START. The two data items associated with this 
command are the time i;ainfall begins on the 
watershed and a code for punching output data. 
If a storm is to be routed through a watershed 
only once, the punch code is deleted. However, 
if more than one routing is to be performed, set 
the punch code equal to a positive number, and 
the output data for the first routing will be 


methods other than the one used in HYMO. The 
input data required for STORE HYD are storage 
location number, hydrograph identification 
number, time increment, watershed area, and 
flow rates of the hydrograph at the specified 
time increment. 

The COMPUTE HYD command is used to 
compute hydrographs from the incremental 
areas of the watershed. The first five items of 
data are storage location number, hydrograph 
identification number, time increment, water- 
shed area, and SCS runoff curve number (.?). 
Normally, data items 6 and 7 are watershed 
height and main stem length. The height and 
length are used to compute the recession con- 
stant K and the time to peak tp . However, if K 
and tp are known or estimated by some other 
method, they can be entered directly into the 
program. This is accomplished by placing a 
minus sign before the values of K and tp and 
entering them as data items 6 and 7, respective- 
ly, The remaining data items are values of the 
mass rainfall at the specified time increment. 

Since most watersheds have a limited number 
of rain gages, the same mass rainfall data may be 
used to develop several hydrographs. Once the 
mass rainfall data have been entered in a 
COMPUTE HYD command, they can be re- 
peated for any number of COMPUTE HYD 
commands without repunching the data. 
Instead, punch a negative number for the eighth 
data item of all COMPUTE HYD commands that 
use the same rain gage. When data from another 
rain gage are entered, the data from the first rain 
gage are lost and cannot be recalled by using the 
negative number code. 

The RECALL HYD command is one of the 
two commands that are computer punched. 
When the punch code is a positive number, the 
output from STORE HYD and COMPUTE HYD 
are punched on c^ds with the RECALL HYD 
command. The RECALL HYD command stores 
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the computed and stored hydrographs on cards; 
it is therefore not necessary to recompute 
hydrographs for future routings. Instead, the 
previously computed hydrographs are read into 
the program, thus saving considerable computer 
time. 

Although the input data for the RECALL 
HYD command are never punched manually, a 
list of the data items may be helpful in checking 
computer-punched cards. The input data are 
storage location number, hydrograph identifica- 
tion number, time increment, drainage area, 
peak flow rate, runoff volume, number of 
hydrograph points, and flow rates of the hydro- 
graph. 

The PRINT HYD command is used to print 
coordinates of a hydrograph, volume of runoff, 
and peak flow rate . The required input data are 
the storage location number and a peak-volume 
code. The peak-volume code is deleted if a 
complete hydrograph printout is desired. If a 
printout of only the runoff volume and the peak 
flow rate is needed, the peak-volume code is set 
to a positive value. 

The PUNCH HYD command is used to punch 
any hydrograph in a HYMO program in the 
proper form for the RECALL HYD command. 
PUNCH HYD has two purposes: (1) If the 
punch code is not used, PUNCH HYD can be 
used to punch one or more hydrographs for 
future use; and (2) if it is desirable to punch 
outflow hydrographs associated with ROUTE, 
ROUTE RESERVOIR, or ADD HYD, PUNCH 
HYD must be used because the punch code only 
provides for punching hydrographs associated 
with STORE HYD and COMPUTE HYD. The 
only datum required for PUNCH HYD is the 
storage location number of the hydrograph to be 
punched. 

The PLOT HYD command is used to plot 
hydro graphs in a HYMO program. It will plot 
one hydrograph on a set of axes, or if a 
comparison is desired, it will plot two hydro- 
graphs on the same set of axes. The required 
input data are the storage location numbers of 
the hydrographs to be plotted. 

The ADD HYD command adds the coordi- 
nates of any two hydrographs. The hydrographs 
are added at a time increment equal to that of 
the hydrograph with the shorter time increment. 
The only data required are the storage location 
number and hydrograph identification number 


of the added hydrograph and the storage loca- 
tion numbers of the two hydrographs to be 
added. 

The STORE RATING CURVE command is 
used to store rating curves that have been 
measured or computed previously. STORE 
RATING CURVE will save considerable com- 
puter time if measured or computed rating 
curves are available. The input data are the 
storage location number, valley, section number, 
and individual rating curve points described by 
elevation, end-area, and flow rate. The number 
of points used to describe a rating curve is 
limited to 20. 

The COMPUTE RATING CURVE command 
is used to compute the stage-area-flow relation- 
ship for a valley section. The input data are 
storage location number, valley section number, 
number of segments in the valley section, 
minimum elevation, maximum elevation, chan- 
nel and flood-plain slopes. Manning’s n value and 
segment boundary point for each segment, and 
horizontal and vertical position of points de- 
scribing the valley section. 

The storage location numbers of the valley 
sections in a particular reach must begin with 1 
and increase by one for each valley section in 
the reach. However, the numbers are assigned 
without regard to upstream or downstream 
order. The valley section identification number 
can be any number from 0.1 to 999.9, These 
rules concerning storage location and valley 
section identification numbers also apply to the 
STORE RATING CURVE command. 

Normally, valley sections are divided into 
three segments (two flood-plain segments and a 
channel segment) for computing the rating 
ciuve. However, some valley sections may have 
more than one channel or may have an extreme 
variation in n values across the flood plain, thus 
requiring more than three segments. A maxi- 
mum of six segments is permitted. Manning’s n 
values for each segment are input with seg- 
ment boundary point (distance from the begin- 
ning of the valley section to the end of the 
segment). Flood-plain n values are positive and 
channel n values are negative. 

Twenty points are used to define a rating 
curve. The location of the points is determined 
by dividing the difference between the maxi- 
mum and minimum elevations into 19 equal 
increments. 


7 



The COMPUTE TRAVEL TIME command is 
used to compute the normal flow travel time 
relationship used in ROUTE. The input data are 
storage location number, reach identification 
number, number of valley sections in the reach, 
reach length, and slope, The reach identification 
number can be any number from 0,1 to 999,9. 
The maximum number of valley sections per 
reach is six. The slope can be either the channel 
or flood-plain slope or a weighted average of the 
two. If flow is confined to the channel, the 
channel slope is of course applicable. If most of 
the flow is in the flood plain, usually the 
flood-plain slope is used. However, a weighted 
slope based on the relative rates of flow in the 
channel and the flood plain may be used. 

The COMPUTE TRAVEL TIME command 
considers each rating curve in the reach in 
computing the travel time flow relationship. 
COMPUTE TRAVEL TIME automatically 
selects the flow rates that are used in computing 
individual travel times. The flow rates of the 
rating curve with the lowest maximum flow rate 
are chosen. If the flow rates of any other rating 
curve in the reach were chosen, the rating curve 
with the lowest, maximum flow rate would have 
to be extrapolated. The travel time table is 
limited to 19 points because of the 20-point 
limit for rating curves. 

The STORE TRAVEL TIME command is one. 
of the two computer-punched commands. When 
the punch code is a positive number, the output 
from COMPUTE TRAVEL TIME is punched on 
cards with the STORE TRAVEL TIME com- 
mand. Therefore, it is not necessary to re- 
compute rating curves or travel time for future 
routings. Instead, STORE TRAVEL TIME reads 
the previously computed travel time flow rela- 
tionship into the program, thus saving consider- 
able computer time. 

The input data for STORE TRAVEL TIME 
are not punched manually, but a list of data 
items may be helpful in checking computer- 
punched cards. The input data are storage 
location number, reach identification number, 


reach length, slope, and individual points of the 
relationship defined by depth, flow, and travel 
time , 

The ROUTE command is used to route floods 
through streams and valleys. The input data are 
storage location number and hydrograph identi- 
fication number of the outflow hydrograph, 
storage location number of the inflow hydro- 
graph, and time increment. The storage location 
number of the outflow hydrograph must be the 
same as the storage location number used in 
COMPUTE TRAVEL TIME for the reach. To 
prevent unnecessary program stoppage, ROUTE 
extrapolates the travel-time table when it is 
exceeded and writes the message, “TRAVEL 
TIME TABLE EXCEEDED.” 

The ROUTE RESERVOIR command is used 
to route floods through reservoirs. The input 
data are storage location number and hydro- 
graph identification number of the outflow 
hydrograph, storage location number of the 
inflow hydrograph, and individual points of the 
reservoir’s outflow-storage relationship. The out- 
flow-storage relationship must be expressed in 
20 points or less. If the outflow-storage relation- 
ship is exceeded, ROUTE RESERVOIR will 
extrapolate the relationship and write the mes- 
sage, “STORAGE-DISCPIARGE TABLE EX- 
CEEDED.” 

The ERROR ANALYSIS command is used to 
determine the error standard deviation and the 
percentage error in peak flow between any two 
hydrographs in a HYMO program. These func- 
tions make ERROR ANALYSIS useful in re- 
search. The input data are the storage location 
numbers of the two hydrographs to be analyzed. 

The SEDIMENT YIELD command is used to 
compute the sediment yield at any point in a 
watershed. Input data required are storage loca- 
tion number of the hydrograph from the area, a 
soils factor, a crop factor, a slope length and 
gradient factor, and a conservation practice 
factor (5). 

The FINISH command is used to end HYMO 
programs. There are no data associated with 
FINISH. 


EXAMPLE HYMO PROGRAM 

A short example problem is presented to Riesel, Tex. A flood will be routed through the 
demonstrate HYMO. Figure 4 is a map of the watershed in its present condition, and the 
6.84-square-mile Brushy Creek watershed near routed outflow hydrograph will be compared to 



the hydrograph measured at gaging station G. 
Also the sediment yield will be predicted and 
compared with the measured sediment yield. 
Then the same flood will be routed through the 
watershed with two proposed reservoirs. To 
determine the effects of the reservoirs, the 
outflow hydrograph and sediment yield will be 
compared to the outflow hydrograph and sedi- 
ment yield of the present-condition routing. 


Comment cards and keywords are used liber- 
ally in the example problem to acquaint the user 
with HYMO. After becoming familiar with 
HYMO, the user may write fewer comments and 
keywords, but generally users find them both 
quite helpful in describing the problem. To save 
space in the example problem, f^ of the 
hydrographs are printed or plotted. Some users 
may choose to print and plot all hydrographs. 


LEGEND 



WATERSHED BOUNDARY 
SU8WATERSHE0 BOUNDARY 
VS- 1 .VALLEY CROSS SECTION NO. I 
1 ROUTING REACH NUMBER I 
UNDERLINED STATIONS NOT IN 
OPERATION 

PROPOSED RESERVOIR SITE 
NUMBER 

INCREMENTAL AREA NUMBER 
RECORDING RAIN GAGE 




SCALE 


I yz 1^4 0 

MILE 


I 

MILE 


Figure 4. Brushy Creek watershed near Riesel, Tex. 
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EXAMPLE HYMO PROGRAM 


a u. 
^ a 
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RUNOFF VOLUME == 4,933 INCHES 

PEAK DISCHARGE RATE = 2360,0 CFS 
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COHPUTe RATING CURVE 10 = 2 VS 7 3 SEGMENTS MIN ELEV=506-'9 FT 

MAX ELEV=515 FT CH SL0PE=.0018 FP SL0PE=-00245 
N==-05 D1ST=636 N=--03 0IST=657 N=.05 0IST=2000 
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BEFORE routing THE TRAVEL TIME - DEPTH - FLOW RELATIONSHIP HOST BE COMPUTED 
FOR THE REACH. 
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THE ROUTE COMMAND 1$ USED TO ROUTE THE HYDRO GRAPH FROM AREA 301 THROUGH REACH* 
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COMPUTE HYD ID=1 HYD N0=302 DT=.16S6667 HR QA=- 837 SQ MI 

CN=82 HEIGHT^69 FT LENGTH=1.33 MI 
REMAINING DATA DEFINE HASS RAINFALL FROM 
START TIME TO END IN INCREMENTS OF OT 0 0 .18 .4 
.89 1.45 2.05 2-7 3-15 3.5 4 4.49 4.72 4.78 4-92 5 


:ki/>p^^^,^^^^^^qqc<occcc>&ocC'OCCocoocci 
o u- 
.J u 
u. 


'4fM<7''0fsl0' ■Or>J0'»0{V0^ir' 


3to<\j^H»-4ooO‘0'0'oo®®r-r''Ovointf>'t-4-fon(<\rv><s^<>ONfM<MM<vjr-< 
o a ^ M «H 
-J o 
UL 



o: 

a 


15 


RUNOFF VOLUME - 4.124 INCHES 

PEAK DISCHARGE RATE = 572.0 CFS 



OUTFLOW HYOROGRAPH REACH 
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IT MAY BE DESIRABLE TO PLOT ALL HYDRQGRAPHS, HOWEVER TO SAVE SPACE IN THIS 
EXAMPLE^ ONLY A FEW OUTFLOW HYORQGRAPHS ARE PLOTTEO- 
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RUNOFF VCLUHE = 4.131 INCHES 

PEAK DISCHARGE RATE = 673.5 CFS 



OUTFLOW HYDROGRAPH REACH 
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RUMOFF VOLUME - 4,196 INCHES 

PEAK DISCHARGE RATE = 4367-5 CFS 



SEDIMENT YIELD ID=6 S0IL=.34 CR0P=-16 CP=.4 LS=.3 
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* THE MEASURED SEDIMENT YIELD FDR THIS FLOOD HAS 3916 TONS* 

* THE ROUTED OUTFLOW HYDROGRAPH CAN BE COMPARED WITH THE MEASURED HYDROGRAPH BY 

* PLOTTING AND ERROR ANALYSIS 
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THE ERROR ANALYSIS COMMAND IS USED TO COMPUTE THE MAGNITUDE OF THE ERROR 
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THE ROUTE RESERVOIR COMMAND IS USED TO ROUTE HYO 301 THROUGH RESERVOIR 502 
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RUNOFF VOLUME ^ 0.382 INCHES 

PEAK DISCHARGE RATE ^ 20. C CPS 



NEXT ROUTE THE OUTFLOW HYD FROM RESERVOIR 501 THROUGH REACH 2. 

BEFORE ROUTING THE TRAVEL TIME - DEPTH - FLOW KELATTONSHIP MUST BE COMPUTED 

for the reach. 
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TO QQTAIN THE OUTFLOW HYD FROM REACH 2 THE HYD FROM AREA 302 MUST &E COMPUTED^ 
AND ADDED TO THE ROUTED HYD 101- 
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OUTFLOW HYOROGRAPH REACH 
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RUNOFF VOLmS = 2,375 INCHES 

PEAK DISCHARGE RATE = 1201.8 CFS 
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ADp THE WD F<^OM AREA 30A TO THE PARTIAL HVD 103 TO OSTAIN THE OUTFLOW FROM 
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OUTFLOW HYOROGRAPH REACH 
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RUNOFF VOLUHH ^ 2*838 INCHES 

PEAK DISCHARGE RATE ^ 2544*1 CFS 



SEDIMENT YIELD 10=6 SCIL=*34 CR0P=,I6 CP=.^ LS=. 3 

SEDIMENT YIELD = 244A-6 TONS 

♦ RECALL THE REACH 5 OUTFLOW HYD ROUTED WITH PRESENT CONDITIONS 
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PLOT THE REACH 5 OUTFLOW HYOS ROUTED WITH PRESENT CONDITIONS AND WITH 
RESERVOIRS- 

PLOT SYMBOLS PRESENT CONDITIONS — G — WITH RESERVOIRS—* — 


o * 
o * 
o + 
o + 
o * 
o + 
o + 


! o + 
•I o + 
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AM»ENDIX 

HYMO MAIN PROGRAM AND SUBROUTINES 


1 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


//SOPTIONS 

C HYMO 

C THIS PROGRAM IS A PROBLEM ORIENTED COMPUTER LANGUAGE FOR 
C DEVELOPING WATERSHED MODELS. 

C 

C 

C DEFINITION OF VARIABLES 

C ID =» STORAGE LOCATION NUMBER 
C NHD * HYDROGRAPH IDENTIFICATION NUMBER 

C DT = TIME INCREMENT IN HOURS 

C DA = DRAINAGE AREA IN SQ Ml 
C OCFS " FLOW RATE ARRAY 

C OCFSF = CUTOFF FLOW RATE IN CFS 
C I END = NUMBER OF POINTS IN A HYOROGRAPH 

C SCFS “ FLOW RATE ARRAY FOR TRAVEL TIME TABLE 

C A = END AREA IN SQ FT 
C Q - FLOW RATE ARRAY FOR RATING CURVES 
C RAIN ^ MASS RAINFALL ARRAY 

COMMON CFS (3001 fOCFSC 300f 6) i IEND(6I fOATA(310) f DA{6 ) t DP( 20 ) ,NPU tNHO 
1, SCFS( 20ltC(20) tA(20,6) tOC 20^6) ♦ RAI NI 2001 t DEEP (20f 6 ) f NER, MAXNO ,NCO 
2MMf ICCf NCOOEf DIST(6) , SEGN(6) , CT8LE ( 50 1 1 1 ) t ITBLE( 50t2) |ZALFA( 20) |DT 
3{ 6)tTIMEfPEAK(6)fROINf I SG(6) 

NCaOE=0 

C READ ZALFA ARRAY 

READ (5r20) ( ZALFA( I ) t 1*1 t 1 5) 

C READ NUMBER OF COMMANDS 

READ (5f21) NCOMM 
C READ COMMAND TABLE 

READ (5,22) ( (CTBLE ( I , J ) , J*1 , 11) , ( I TBtl: ( I , J ) , J* 1, 2 ) , I = 1,NC0MM ) 
WRITE (6,23) (ZALFA(n,I = l,15) 

WRITE (6,24) 

WRITE (6,2 5) ((CTBLE(I,J) ,J=:1,11),(ITBLE(! , J ) , J*l , 2 ) , I *1 , NCOMM ) 

C ZERO CODES — NPU * PUNCH, ICC * CONTINUATION CARD, NER = ERROR 
NPU*0 
ICC*0 

1 NER*0 
CALL HONDO 

IF (NER) 2,2,19 

2 GO TO (3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19), NCODE 

3 TIME^DATAd) 

NPU=:0ATA(2) 

GO TO 1 

4 CALL STHYD 

GO TO I 

5 CALL RECHD 

GO TO 1 

6 CALL CMPHYO 

GO TO 1 

7 CALL PRTHYD 

GO TO 1 

8 CALL PUHYO 

GO TO 1 

9 CALL HPLOT 

GO TO I 

10 CALL ADHYO 

GO TO I 

11 CALL SRC 
GO TO 1 

12 CALL CMPRC 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
IT 
la 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 
4B 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


54 



35 


GO TO 1 

A 

60 

36 

13 

CALL STT 

A 

61 

37 


GO TO 1 

A 

62 

30 

14 

CALL 'CMPTT 

A 

63 

39 


GO TO 1 

A 

64 

^0 

15 

CALL ROUTE 

A 

65 



GO TO 1. 

A 

66 

^2 

16 

CALL RESVO 

A 

67 

43 


GO TO 1 

A 

66 

44 

17 

CALL ERROR 

A 

69 

45 


GO TO 1 

A 

70 

46 

18 

CALL SEDT 

A 

71 

47 


GO TO 1 

A 

72 

48 

19 

STOP 

A 

73 


C 


A 

74 

49 

20 

FORMAT «15AIJ 

A 

75 

50 

21 

FORMAT (12) 

A 

76 

51 

22 

FORMAT (2Al,9A2t2I3) 

A 

77 

52 

23 

FORMAT ( iHl,9Xt8HZALFA = ,15AI///) 

A 

70 

53 

24 

FORMAT (16X,13HC0HMAND TABLE//) 

A 

79 

54 

25 

FORMAT ( 10X»2AI ,9A2,2I3) 

A 

80 

55 


END 

A 

81- 

56 


SUBROUTINE HONDO 

8 

1 


C 

THIS SUBROUTINE READS IN A DATA CARO , SEARCHES AN ALPHAMERIC 

0 

2 


C 

CODE TABLE TO DETERMINE THE NCODE OF THE OPERATION AND 

B 

3 


c 

COLLECTS VARIABLES FROM THE FREEFLOATING DATA FIELD 

B 

4 

5.7 


COMMON CFS (300) tOCFS( 300,6) , lENO (6 ) ,OATA( 310) , 0A(6 ) , 0P( 20 ) , NPU ,NHD 

B 

5 



l,SCFS(20) ,C(20) , A(20,6) ,0(20,6) ,ftAIN (200) ,DEEP( 20,6) ,NER,MAXNO,NCO 

B 

6 



2MM,ICC,NCODE,OIST(6),SEGN(6),CTBLE{50,ll),ITBLe(50,2) ,ZALFA ( 20 ) , DT 

B 

7 



3(6),TIME,PEAK ( 6 ) , RO IN , 1 SG ( 6 ) 

3 

8 

50 


DIMENSION CHAR(60), ALPHA(ll) 

B 

9 

59 


DIMENSION AUXA(IO), AUXBIIO) 

B 

10 

60 


IF (ICC) 1,1,3 

a 

11 


c 

READ IN DATA CARO 

B 

12 

61 

1 

READ (5,A2) ( ALPH A( 1) , I = l , 1 1 ) , ( CHAR( I ) , I = l ,60 ) 

B 

13 


c 

IF FIRST CHARACTER IS BLANK THE CARD IS A CONTINUATION OF 

B 

14 


c 

PREVIOUS CARD. 

B 

15 

62 


IF (ALPHA! 1)-ZALFA( ID) 2,9,2 

B 

16 

63 

2 

IF ( ICC) 3,3,40 

B 

17 


C 

ASTERISK IN COL. 80 MEANS SKIP TO NEW PAGE BEFORE PRINTING CARO 

B 

18 

64 

3 

IF (CHAR(60)-ZALFA( U) ) 4,5,4 

B 

19 

65 

4 

WRITE (6,43) 

B 

20 

66 

5 

WRITE (6,44) (ALPHA(!),I«1,1I),{CHAR(I),I=1,60) 

B 

21 


C 

IF FIRST CHARACTER IS A * THE PREVIOUS CARD WAS A COMMENT CARD 

& 

22 

67 


IF (ALPHA(1)-ZALFA(12)) 10,6,10 

B 

23 


C 

IF PUNCH CODE POSITIVE, COMMENT CARDS ARE PUNCHED. 

B 

24 

68 

6 

IF (NPU) 8,8,7 

B 

25 

69 

7 

WRITE (7,45) (ALPHA(I),I=l,ll),(CHAR(I),l=l,60) 

B 

26 

70 

8 

ICC=0 

B 

27 

71 


GO TO 1 

6 

20 

72 

9 

WRITE (6,44) (ALPHA(I),I=l,ll),(CHAR(l),I«l,60) 

B 

29 

73 


GO TO 24 

B 

30 


C 

SEARCH FIRST TWO ALPHAMERIC CHARACTERS TO SEE IF THEY ARE NUMBERS 

B 

31 

74 

10 

ICC = l 

B 

32 

75 


DO 12 1=1,10 

B 

33 

76 


IF (ALPHA! l)-2ALFA(I)) 11,15,11 

B 

34 

77 

U 

IF ( ALPHA! 2)-ZALFA( I) ) 12,15,12 

B 

35 

78 

12 

CONTINUE 

B 

36 


C 

STATEMENT NUMBER 7 IS BRANCHED TO IF NUMBERS ARE PRESENT 

B 

37 
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c 

c 


79 


80 


81 

C 

82 

13 


C 


C 

83 


84 


85 


86 

14 


C 


C 

87 


88 


89 

C 

90 

15 


C 

91 


92 


93 

16 

94 


95 

17 


C 


C 

96 


97 


98 

18 

99 


100 

19 

101 

20 

10 2 

21 

10 3 

22 


C 

104 

23 

105 

47 

106 


107 

24 

108 

25 

109 


no 

26 

111 


112 


113 


114 


115 

C 

116 

27 


C 

117 


118 

28 

119 


120 

29 

121 


122 

30 


C 

123 

31 


[F NOT NUMBER SEARCH COMMAND TABLE FOR MATCH 
CALL FIRST 10 VALUES FROM PERMANENT DATA STORAGE 
DO 14 I=lf50 
00 13 J=l,ll 

IF (CTBLE( I,J)-ALPHA(J) ) I4,l3fl4 
SN lO^PART MATCH 

CONTINUE ^ 

IF THIS LOOP IS COMPLETED WE HAVE COMPLETE MATCH- CALL NCOOE 

AND MAX NUMBER AND EXIT LOOP 

NCODe=IT0LE(I 

MAXNO^ITBLEI 1 #2) 

GO TO 21 
CONTINUE 

IF MAJOR LOOPS FINISHED WITHOUT A MATCH WRITE ERROR MESSAGE 

AND SET NER = 1 

NER=:1 

WRITE <6,46) 

RETURN 

CONVERT DIGIT INPUT CODE FROM ALPHAMERIC TO INTEGER FORM 
NCODE=GIT< ALPHA 

FIND MAX NUMBER OF DATA ITEMS FOR THIS NCOOE 
DO 17 Ul,50 

IF ( IT6LE( I ,n-NCOOE) 17,16|17 
MAXN0^ITBLE(1#2) 

GO TO 21 
CONTINUE 

SEARCH DATA ROUTINE 

SEE IF ANY DATA FOR THIS CARD 

DO 19 I«l,50 

IF (ITBL6( I,1)-NCDDE) 19,18,19 
MAXNO=ITBLE(I ,2) 

GO TO 20 
CONTINUE 
CONTINUE 

IF (MAXNO) 23,22,23 
RETURN 

ZERO ARRAYS AND COUNTERS 
DO 47 1-1,310 
DATA ( 

NDATA=1 
NCHAR-0 
DO 26 1^1,10 
AUXA(I)^0, 

AUXB( I)=0. 

IT 1=1 
IT2=1 
SIGN=l. 

L0GIT=0 

KDGIT=0 

CARRY OUT DIGIT BY DIGIT SEARCH AND ACCUMULATION 
NCHAR=NCHAR+1 

HAVE WE CONSIDERED ALL CHARACTERS - RETURN IF SO 
IF (NCHAR-60) 28,32,1 
DO 29 1=1,15 

IF (CHARINCHARI-ZALFAID) 29,30,29 

CONTINUE 

GO TO 32 

GO TO <33,33,33,33,33,33,33,33,33,33,32,27,36,32,31,27), I 

SN 39 HANDLES SIGN CONTROL ON 1130 VERSION 

SIGN=-KO 


6 38 

B 39 
B 40 
B 41 
B 42 
B 43 
B 44 
B 45 
B 46 
B 47 
B 48 
B 49 
B 50 
B 51 
B 52 
B 53 
B 54 
0 55 

0 56 

B 57 
0 50 

0 59 

8 60 
B 61 
B 62 
B 63 
D 64 
0 65 

B 66 
B 67 
B 68 
B 69 
B 70 
B 71 
B 72 
B 73 
B 74 
0 75 

B 76 
B 77 
B 78 
B 79 
B 80 
B 81 
B 82 
B 83 
B 84 
B 85 
B 86 
B 87 
B 88 
B 89 
B 90 
B 91 
B 92 
B 93 
B 94 
B 95 
B 96 
B 97 
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12^ 


GO TO 27 

6 

98 


C 

CHARACTER IS BLANK OR COMMA - DOES IT FOLLOW A DIGIT 

B 

99 

125 

32 

GO TO (27,48) , ITl 

B 

100 


C 

character is a digit - HAS A DECIMAL BEEN ENCOUNTERED 

B 

101 

126 

33 

GO TO ( 34,35) , IT2 

B 

102 

127 

34 

LDGIT=LDGI T+1 

B 

103 

128 


1 Tl=2 

B 

104 

129 


AUXAILDGIT )=:CHAR(NCHAR) 

6 

105 

ISO 


GO TO 27 

8 

106 

131 

35 

K0GIT=K0GIT4-1 

B 

107 

132 


AUXB(KDG1T )==CHA(UNCHAR) 

B 

108 

133 


GO TO 27 

B 

109 


C 

CHARACTER IS A DECIMAL - DOES IT FOLLOW A DIGIT 

B 

no 

134 

36 

GO TO ( 37,38) , ITl 

B 

111 

135 

37 

I Tl-2 

8 

112 

136 


LOGIT^l 

B 

113 

137 

38 

IT 2 = 2 

8 

114 

138 


GO TO 27 

B 

U5 


C 

ROUTINE TO CONVERT ALPHABETIC ARRAY TO FLOATING POINT NUMBER 

B 

116 

139 

48 

DATA (NDATA)^GIT(AUXA,1,LDGIT,K)^GIT(AUX0, 1,10,0.) 

8 

117 

140 


DATA ( NDATA)=:DATA(NDATA)*SIGN 

B 

118 


C 

IS ALL DATA FURNISHED YES-RETURN NO INCREASE N DATA KEEP ON 

B 

119 

141 


IF ( NDATA-MAXNO) 41,39,39 

6 

120 

142 

39 

ICC^O 

B 

121 

143 

40 

RETURN 

B 

122 

144 

41 

NDATA=NDATA+l 

B 

123 

145 

C 

GO TO 25 

B 

B 

124 

125 

146 

42 

FORMAT (2A1,9A2,60A1) 

B 

126 

147 

43 

FORMAT (IHl) 

•B 

127 

140 

44 

FORMAT (5X ,2Al,9A2,60Al) 

B 

128 

149 

45 

FORMAT (2A1 ,9A2,60Al) 

B 

129 

150 

46 

FORMAT (IOX,20HCOMMANO NOT IN TABLE) 

B 

130 

151 


END 

B 

131- 

152 


FUNCTION GIT ( TCARO , J , JLAST, SHI FT) 

C 

1 

153 


DIMENSION TCARD(IO), A(10) 

C 

2 

154 


DATA A( 1)/1HI/,A(2)/IH2/,A(3)/IH3/, A(4) /IH4/, A< 5) /IH5/, A(6) /H 6/ 

c 

3 

155 


DATA A(7)/IH7/,A(8) /I H8/ , A( 9) /IH9/ , A( 10) /IHO/ 

c 

4 

156 


GIT^O. 

c 

5 

157 


TEN=10. 

c 

6 

158 


SUM=0. 

c 

7 

159 


DO 3 JNOW=J, JLAST 

c 

8 

160 


TTEST=TCARD( JNOW) 

c 

9 


C 

CHECK FOR LAST ENTRY 

c 

10 

161 


IF (TTESr.EQ.O.) GO TO 4 

c 

n 


C 

FIND NUMBER AND COMPUTE VALUE 

c 

12 

162 


DO 2 NUMB=l,lO 

c 

13 

163 


IF (TTEST-A(NUMB)) 2,1,2 

c 

14 

164 

1 

ZTEST=NUM0 

c 

15 

165 


IF (ZTEST.EQ.IO.) ZTEST^O. 

c 

16 

166 


SUM==SUM*TEN+ZTEST 

c 

17 

167 


GO TO 3 

c 

18 

163 

2 

CONTINUE 

c 

19 

169 

3 

CONTINUE 

c 

20 

170 

4 

IF (SHIFT) 6,5,6 

c 

21 

L7l 

5 

FI^JNOW-1 

c 

22 

172 


SUM=SUM*(0.1>^*FI) 

c 

23 

173 

6 

GIT=SUM 

c 

24 

174 


RETURN 

c 

25 
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175 


26 - 




176 


SUBROUTINE STHYD 

D 

1 


C 

THIS SUBROUTINE STORES THE COORDINATES OF HYOROGRAPHS. 

0 

2 

177 


COMMON CFS I 300) tOCFS ( 300 1 6) , lENOl 6 ) f DATA! 3 10) t DA 1 6 ) , DP ( 20 ) t NPU,NHO 

D 

3 



li SCFSI 20 ) »C( 20) ,A( 20»6) tO (20*6) i RAIN (200 ) , DEEP! 20* 6 ) * NER * MA XNO * NCO 

D 

4 



2HM ,ICCf NCQDE , 0IST(6 ) *SEGN<6) ,CTBLEl50*ll)* ITBLE(50*2) *2ALFA(20)*DT 

D 

5 



3( 6) ,TIME,PEAK 16) *R0IN*!SG(6) 

0 

6 

178 


ID=DATA(1) 

D 

7 

179 


NHD=DATA(2) 

0 

8 

180 


0T(ID)=0ATAI3) 

D 

9 

181 


DA(ID)=DATA<4) 

D 

10 

182 


J = 5 

D 

11 


C 

REMAINING DATA ARE FLOW RATES 

D 

12 

183 


OCFSIl , ID) =OATAU) 

0 

13 

18A 


PEAKd 0) = 1 . 

D 

14 

185 


RO=DATA( J) 

D 

1 5 

186 


DO <t 1=2,300 

D 

16 

187 


J = J+1. 

D 

17 

188 


OCFS( I , ID) =OATA( Jl 

0 

18 

189 


RO=RO+OCFS(1, ID) 

0 

19 


C 

IS FLOW RECEDING 

0 

20 

190 


IF (OCFSn,IO)-OCFSII-l*IO)» 1,2,2 

0 

21 


c 

HAS FLOW RECEOEO TO CUTOFF RATE 

D 

22 

191 

1 

IF (OCFSIl, IDI) 5,5,6 

D 

23 


c 

DETERMINE PEAK FLOW 

D 

24 

192 

2 

IF (OCFSIl ,ID)-PEAK(IOI ) 6,4,3 

D 

25 

193 

3 

PEAKI!OI=OCFS( I ,10) 

D 

26 

19A 

4 

CONTINUE 

0 

27 

195 

5 

IEND(1DI=1-1 

D 

28 

196 


M=IENO( ID) 

D 

29 

197 


ROIN’I R04DT( ID) )/IDA( 101*665,333) 

D 

30 


C 

PUNCH CODE 

D 

31 

198 


IF (NPU) 7,7,6 

D 

32 

199 

6 

WRITE (7,8) ID,NHD,0T(1D),0A(I0),PEAK(1D),R0IN,IEND(ID> 

D 

33 

200 


WRITE (7,9) (OCFS(J,IO) ,J=1,M) 

D 

34 

201 

7 

RETURN 

D 

35 


C 


0 

36 

202 

6 

FORMAT! ‘RECALL HYO • , T2 1 , ' I 0= • , U , T29 , • HYD N0= • , I 3 , T62 , ' 0T= • , F9. 

D 

37 



16,* HRS' ,T61 , 'DA=* ,F8.3,* SO MI • /T2 1 , ' PE AK= ' , F7. 0 , ' CFS ' ,T60 , • R0=' , 

D 

38 



2F6.3,' INCHES', T59, 'NO PTS= ' , I 3/T2 1 , ' FLOW RATES') 

D 

39 

203 

9 

FORMAT (T2l,7F0.O> 

0 

40 

204 


END 

D 

41- 

20 5 


SUBROUTINE RECHD 

E 

1 


C 

THIS SUBROUTINE RECALLS PREVIOUSLY COMPUTED AND PUNCHED 

E 

2 


C 

HYDROGRAPHS 

E 

3 

206 


COMMON CFS (300) ,OCFS( 300, 6) , I ENO( 6) ,OATA( 3 10 ) , DA( 6 ) ,DP( 20) ,NPU,NHD 

E 

4 



1, SCFSI 20) ,C(20) , AI 20, 6) ,0(20,6) , RAIN! 200) , DEEP ( 20 , 6 ) , NER , MAXNO, NCO 

E 

5 



2MM,iCC,NCDDe,0IST(6) ,SEGN(6),CT0LE(50,ll) , ITBL E{ 50 , 2 ) , ZALF A ( 20 ) , DT 

E 

6 



3( 6),T1ME,PEAK(6) ,R01N, 1S6(6) 

E 

7 

207 


1D=0ATA(1) 

E 

8 

208 


NH0=DATA(2 ) 

E 

9 

209 


DT(I0)=DATA(3) 

E 

10 

210 


0A(1D)=DATA(6) 

E 

11 

211 


PeAK(ID)=DATA(5) 

E 

12 

212 


R0IN=DATAI6) 

E 

13 

213 


IEND(ID)=0ATA(7) 

E 

14 

214 


M=IENO(IO) 

E 

15 

215 


J = 0 

E 

16 




C 

REMAINING DATA ARE FLOW RATES 

E 

17 

216 


00 1 I=1»M 

E 

18 

217 


OCFS( I , ID)=DATA( Jl 

E 

19 

218 

1 

J = J+1 

E 

20 

219 


RETURN 

E 

21 

220 


END 

e 

2 2- 

221 


SUGftOUTINE CHPHYO 

F 

l 


C 

THIS PROGRAM DEVELOPS A UNIT HYDROGRAPHt CONVERTS MASS RAINFALL 

F 

2 


c 

TO POINT RUNOFFf AND COMPUTES STORM HYDROGRAPHS BY SUMMATION. 

F 

3 

222 


COMMON CFS (300) ,0CF5( 300,6) 1 1 END ( 6 ) , DATA ( 3 10 ) , 0A( 6 ) , DP ( 20 ) ,NPU,NHD 

F 

4 



1, SCFSf 20),C(20),A(20,6)fQ(20,6) ,RAIN(200) , DEEP (20,6) , NER » MAXNO i NCO 

F 

5 



2MM,ICC,NC0DE, DIST(6) ,SEGN(6) tCTBLE (50,11) , I TBL E ( 50 , 2 ) , ZALFA ( 20 ) , OT 

F 

6 



3( 6) ,TIME,PeAK(6) ,ROIN, ISG(6) 

F 

7 

223 


tD=DATA( 1) 

F 

8 

224 


NH0=0ATA(2 ) 

F 

9 

225 


DT( I0)^DATA(3) 

F 

10 

226 


DA( I0)=0ATA(4) 

F 

1 L 

227 


CN-DATA(5) 

F 

12 


c 

ARE K AND TP FURNISHED OR WILL THEY BE COMPUTED 

F 

13 

228 


IF (DATA(6)) 1,2,2 

F 

14 

229 

1 

XK^-DATA(6) 

F 

15 

230 


TP=:-0ATA(7) 

F 

16 

231 


GO TO 3 

F 

17 

232 

2 

HT=0ATA{6) 

F 

IB 

233 


XL=DATA{7) 

F 

19 

234 


SLOPE=HT/XL 

F 

20 

235 


XLOW=( XL«*2.) /DA( 10) 

F 

21 

236 


XK = 27.0*( OA(IO)««.231)«(SLOPE*’«‘(-.7 77) )*(XLDW**.124) 

F • 

22 

237 


TP=4#6 3«( 0A.( I 0) 422 )»)«( SLOPED* (-.46) ) * I XL DW*« . 133 ) 

f 

23 

238 

3 

PEAK! ID)=^1. 

F 

24 

239 


DO 4 1^=1,300 

F 

25 

240 

4 

QCFSI I f ID)==0. 

F 

26 


C 

COMPUTE N BY ITERATION. 

F 

27 

241 


XN=^5.0 

F 

28 

242 


XKTP^XK/TP 

F 

29 

243 


DO 6 1=1,50 

F 

30 

244 


TINF=l.+SQRT( l./(XN-l.) ) 

F 

31 

245 


XNl=.05/( XKTP’frI ALOG( TINF/(TlNF + .05) ) + .05))+l. 

f 

32 

246 


DIFF=AI3S(XM-XN} 

F 

33 

247 


IF (DIFF-.OOl) 7,7,5 

F 

34 

248 

5 

XN=XN1 

F 

35 

249 

6 

CONTINUE 

F 

36 

250 


WRITE (6,29) 

F 

37 

251 


GO TO 28 

F 

38 


C 

DETERMINE Cl. 

F 

39 

252 

7 

DELT=TfNF/100. 

F 

40 

253 


TC1=0. 

F 

41 

254 


XN1P=XN-1. 

F 

42 

255 


XNIM=1.-XN 

F 

4 3 

256 


DO 8 1=2,101 

F 

44 

257 


TCI=TC1+D6LT 

r 

45 

258 

8 

CFS( 1) =(TCL’«^«XNIP)«EXP(XNIM*(TCI-1.)) 

F 

4 6 

259 


SUM=CFS( 10 l)/2. 

F 

47 

260 


DO 9 1=2,100 

F 

48 

261 

9 

SUM = SUM<*CFS( I ) 

F 

49 

262 


C1=$UM*DELT 

F 

50 

263 


CFSI I = CFS( 101) 

F 

51 

264 


TT INF^TINF^^'TP 

F 

52 

265 


TRECl=TTINF+2.«XK 

F 

53 


59 



266 


267 


268 

C 

2 69 


270 


27L 


2 72 

C 

2 73 


274 


275 


276 


277 

10 

278 


2 79 


2 80 


281 


282 

11 

233 

L2 

284 

13 

2 85 


286 

33 

287 


288 


289 

14 

290 


291 


292 

15 


C 

293 


294 


295 


296 


297 


298 

16 

299 


300 

17 

301 

18 

302 


303 

19 

304 


30 5 

20 

306 


307 

21 


C 

308 


309 


310 


311 

22 

312 

23 

313 


314 


315 


316 

24 

3X7 


318 


319 

C 

320 



Eee=EXP((TTINF-TRECi»/XKI 

XK1=3.#XK 

B=6<i5.333/(C1+CFSII*( XKTP*( l.-EEe)+EEE*C XKl/TP > n 
COMPUTE El, QP, AND CFSI. 

QP*(B*0A( 10) »/TP 
CF5!=QP*CFS(101) 

CF SRl=CFSI*EEE 

WRITE (6,30) XN,QP 

DETERMINE INCREMENTAL RUNOFF. 

R= lOOO./CN-lO. 

B1=.2*R 

J*8 

IF (OATA(J)) 13,10,10 
RArN(l)«OATA( J) 

DO 11 1=2,300 
J = J+1 

RA IN(I )=DATA( J) 

IF IRAIN( I )-RAIN( I-U) 12,11,11 

CONTINUE 

NUNB=1-1 

DO 15 1=1, NUMB 

IF (RAIN(I)-Bl) 33,33,14 

DATA ( 1)=0. 

Ql=0. 

GO TO 15 

Q2=(IRAIN( I )-B1)**2.)/(BAIN(I)+.8*R) 

DATA I n = Q2-Ql 

Q1»Q2 

CONTINUE 

COMPUTE UNIT HYDROGRAPH. 

T2=0. 

CFSI 1)=0. 

DO 20 1=2,300 
T2=T2tOTI ID) 

IF (T2-TTTNF) 16,16,17 

CF SI I)=QP*I IT2/TP)**XNIP1*EXP( XNIM* I T2/TP- 1 . I ) 

GO TO 20 

IF IT2-TREC1) 18,18,19 
CFS(1)=CFS1*EXP((TTINF-T2)/XK) 

GO TO 20 

CFSI I)=CFSR1*EXPI ITREC1-T21/XK1I 
IF ICFS(I)-1.) 21,21,20 
CONTINUE 
1 = 300 
ICND=I 

COMPUTE STORM. HYDROGRAPH. 

DO 24 J=2,NUM8 
N=J+lCND-2 
IF IN-300) 23,23,22 
N=300 
KK=J 
1=2 

00 24 K=KK,N 

OCFS(K,lD)=OCFSlK,lD)+DATAIJ)*CFSn ) 

1 = 1 + 1 
M=K-l 
R0=0. 

DO 26 1=2, M 
COMPUTE RUNOFF VOLUME 
RO = RO+OCFSII, ID) 


F 

54 

F 

55 

F 

56 

F 

57 

F 

58 

F 

59 

F 

60 

F 

61 

F 

62 

F 

63 

F 

64 

F 

65 

F 

66 

F 

67 

F 

68 

F 

69 

F 

70 

F 

71 

F 

72 

F 

73 

F 

74 

F 

75 

F 

76 

F 

77 

F 

78 

F 

79 

F 

80 

F 

81 

F 

82 

F 

83 

F 

84 

F 

85 

F 

86 

F 

87 

F 

88 

F 

89 

F 

90 

F 

91 

F 

92 

F 

93 

F 

94 

F 

95 

F 

96 

F 

97 

F 

98 

F 

99 

F 

100 

F 

101 

F 

102 

F 

103 

F 

104 

F 

105 

F 

106. 

F 

107 

F 

108 

F 

109 

F 

110 

F 

111 

F 

112 

F 

113 


60 



c 

321 

3e2 2 5 

323 26 

324 

325 

326 
32 7 

C 

328 

329 27 

330 

331 28 
C 

332 29 

333 30 

334 31 


335 32 

336 

33 7 

C 

33 8 


339 

340 

C 

341 

342 1 

343 

344 2 

345 3 

346 

347 4 

348 5 

349 

350 6 
C 

351 7 

352 e 

353 

354 

355 

C 

356 

357 

358 9 

359 

360 

361 10 

362 

363 

364 11 

365 

366 12 


DETERMINE PEAK 

IF (OCFSU , rDl-PEAKlIon 26,26,25 
PEAKtIDI=QCFS< I ,ID) 

CONTINUE 

I=M 

I END! rOI^I 
M= I 

ROlN^CRO^i'DKIOll/COAt ID)» 64 5^3331 
PUNCH CODE 
IF (NPU) 28,28,27 

WRITE {7,31) 10,NHD,DTnD),DAnD>,PEAKlID> ,ROI N, I END( IDJ 
WRITE <7,32) (aCFS( I, 10) 

RETURN 

FORMAT (*N DID NOT CONVERGE AFTER 50 ITERATIONS.*) 

F ORM AT ( TIO , •SHAPE CONSTANT, N = • , F 6 . 3/T 10 , • UN I T PEAK *,F10.l,*C 
IFSV) 

FORMAK •RECALL H YO » , T2 1 , * 1 D- • , I 1 ,T29, • HYO NO-* , I 3 ,T42, •DT - • ,F9. 
16,* HRS*tT61, •DA^* , F8.3, • SQ MI * /T21 , ' PEAK- * , F7, 0 , • CPS* , T40 , •RO-* , 
2F6.3,* INCHES* ,T59, ‘NO PTS“ * , 1 3/T21 , * FLOW RATES*) 

FORMAT <T21,7F8.0> 

END 


F 114 
F 115 
F 116 
F 117 

F ue 

F 119 
F 120 
F 121 
F 122 
F 123 
F 124 
F 125 
F 126 
F 127 
F 12B 
F 129 

F 130 
F 131 
F L32 
F 133 
F 134” 


SUfiROUTlNE PRTHYO 

THIS SUBROUTINE PRI NTS- THE COORD INA TES OF A HYDROGRAPH* 

COMMON CFS{300),OCFS{300,6),IENO(6),OATA(310),DA<6),DP(20) ,NPU,NHO 
L,SCFS<2D),C<203,A(20,6),Q(20,6), RAIN (200) , DEEP ( 20 • 6 ) , NER , MAXNO , NCO 
2MM ,I CC iNCOOE ,D 1ST (6 ) , SEGN(6 ) ,CTBLE ( 50 , 1 1 1 , ITBLE ( 50 , 2 ) ,ZALF A ( 20 ) , OT 
3( 61 ,TIME,PeAK(6) ,R0IN,1SG<6) 

IO=DATA( 11 
NPK-DATA(2 ) 

DETERMINE TYPE OF HYOROGRAPH 
IF (NHO-100) 6,6,2 


WRITE (6,14) 
GO TO 7 
IF (NHD”3aO) 
WRITE (6,15) 
GO TO 7 
IF (NHD-500) 
WRITE (6,16) 
GO TO 7 
WR ITE <6,17 ) 


NHD 

3.3.4 
NHD 

1.1.5 

NHD 

NHD 


POSITIVE NPK MEANS PRINT ONLY PEAK AND VOLUME 
IF (NPK) 0,8,11 

WR ITE (6,18) 

IEND( ID) 

TI MEUTIME 

BUILD TIME ARRAY IN DATA 
DO 9 1=^1, M 
DATA ( 1) = TIME I 
TIME1=TIME1+DT( ID) 

M4=M+4 

M5=M4/5 


J^J+1 

WRITE (6,19) (0ATA( U ,QCFS( I ,ID) ,I = J,M,M5) 

IF (J”M5) 10,11,11 

WRITE (6,20) RaiN,PEAK<ID) 

IF (NPU) 13,13,12 
WRITE (7,21) ID, NPK 


G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

6 

G 

G 

G 

G 

G 

G 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 


61 



367 13 
C 

368 14 

369 15 

370 16 

371 17 

372 18 


373 19 

374 20 

375 21 

376 


RETURN 


FORMAT 
FORMAT 
FORMAT 
FORM AT 
=ORMAT 
riHE 
IX,U3HHRS 
CFS 


, 13/) 
4/1 


(1H0|46X,21HHYDRGGRAPH FROM AREA 
{lH0,41X,19HPARtlAL HYDROGRAPH ,14/1 
( 1H0,39X,29HQUTFLGW HYDROGRAPH RESERVQ 
aH0,44X,25HQUTFL0W HYDROGRAPH REACH , 
(10X,114HTIME Flow TIME 


CFS 

HRS 


HKb CFS 

IRMAT {51 5X,F10,3,F10.0) ) 

RMAT ( 1HO,9X,16HRUNOFF VOLUME 
ARGE RATE = iF10.1.4H CFS///) 


HRS 


TIME 

CFS 

HRS 


IR I 

14/1 


RATE 

)RMAT( tDOTMT 


I AonrMji'iurr vuuyric 

,F10.1,4H CFS///) 

MT LlV/nt TOl tTn.., 


F10.3,8H INCHES 


14/1 

FLOW 

TIME FLOW/1 

HRS 

CFSl 

/10X,22HPEAK DIS 


G 36 
G 39 
G 40 
G 41 
G 42 
G 43 
G 44 
G 45 
G 46 
G 47 
G 48 
G 49 
G 50 
G 51 
G 52- 


377 

C 

C 

378 


379 

380 

381 

382 

383 

C 

384 1 


385 2 

386 


SUBROUTINE PUHYD H I 

THIS SUBROUTINE PUNCHES HYDROGRAPHS IN FORM TO BE USED BY H 2 

SUBROUTINE RECHD H 3 

common CFS (300) ,OCF S( 300 , 6 ) , I END ( 6) , DAT A ( 3 10 ) » DA ( 6 ) , DP ( 20 ) , NPU , NHD H 4 
1,SCFS( 201,C(20) ,A(20,6) ,Q (20,6 ) , RA! N(200 ) i DEEP ( 20 , 6 ) , NER, MAXN0,NCD H 5 
2MM,ICC,NC00E,DlST(6),SeGN(6) ,CTBLE (50 , 1 1) , I TBlE ( 50, 2 ) , ZALFA ( 20 ) , DT H 6 
3(6), TIME, PeAK(6),R0IN,ISG(6) H 7 

lO^DATA(l) H 8 

M^IENDdO) H 9 

WRITE (7,1) IO,NHO,DT(IO),DA{ID),PEAK(1D),ROIN,IEND(ID) H 10 

WRITE (7,2) (OCFSdtID) ,I = 1,M) H 11 

RETURN H 12 

H 13 

FaRMAT( ^RECALL HYD • , T2 1 , M 0= • , U , T29 , • H YD NO^ • , I 3 , T42 , • DT= • , . H 14 

16, » HRSST61 , *0A==dF8.3,* SQ MP /T21 , ' PEAK= • , F7.0 , ‘CFS • ,T40 , » RO^* , H 15 
2F6.3,* 1NCHESST59, ‘NO PTS= » , 1 3/T2 1 , • FLOW RATESM H 16 

FORMAT {T21,7F8a0) H 17 

END H 18- 


387 

C 

388 


389 

390 

391 

392 

393 

394 

395 

C 

396 

C 

397 1 

398 

399 2 

400 3 

401 

402 4 
C 
C 

403 5 


SUBROUTINE HPLOT 1 

THIS SUBROUTINE PLOTS EITHER 1 OR 2 HYDROGRAPHS ON A SET OF AXIS I 

COMMON CFS (300) ,OCFS( 300,6) , I END ( 6 ) , DATA ( 3 10 ) , DA ( 6 ) ,DP( 20 ) ,NPU ,NHD I 

1, 5CFS( 20) ,C(20) , A( 20, 61 ,0(20,6) , RAIN (200) , DEEP (20, 6) , NER , MAXNO ,NCO I 

2MM,ICC,NC0DE,0IST(6) ,SEGN(6) tCTBLE (50, 11), I TBLE (50,2) , ZALFA (20) ,DT I 

3( 6),TIME,PEAK(6) ,R0IN,ISG(6) 1 

lOl^DATA(l) I 

I02-DATA(2) I 

DATA ZERO, PLUS, BLANK, DASH, DOT/ *0 • , * + ‘ ‘.V I 

MRTO-1 I 

XMRTO=U I 

MAX=U8 T 


ARE THERE i OR 2 HYDROGRAPHS 
IF dD2) 1,1,2 

DETERMINE HIGHEST PEAK IF 2 HYDROGRAPHS 
QHAX=:PEAK( IDl) 

GO TO 14 

IF (PEAK(lDU-PEAK(ID2n 3,3,4 
QMAX=PEAK( ID2) 

GO TO 5 

QHAX=PEAK( 101) 

IF 2 HYDROGRAPHS DETERMINE LARGEST DT AND INTERPOLATE OTHER 
HYDROGRAPH IF NECESSARY 
IF (DT( 101 )-DT( 102) ) 6,13,7 


1 

2 

3 

4 

5 

6 
7 
a 

9 

10 

IL 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


62 



404 

6 

40 5 


406 


407 

7 

408 


409 

8 

410 


411 


412 


413 


414 


415 

9 

416 


417 


418 


419 

10 

420 


421 


42 2 

11 

42 3 


424 


42 5 


426 

12 

427 

13 

42 8 

14 

42 9 


430 

15 

431 

16 


C 

432 

17 

433 


434 


435 

18 

436 

19 


C 

43 7 


43 0 

20 

439 


440 


441 


442 


443 


444 

21 

445 

22 

446 


44 7 

23 

448 

24 

449 

25 

450 


451 


452 


453 


454 

26 

455 

27 

456 

28 

457 


458 


459 

29 

460 



L=I01 
K= ID2 
GO TO 8 
L= ID2' 

K=ID1 
M-IEND(LI 
TID=^DT« Kl 
T I DH=0 . 

DO ll 

TIDH^TIDH+DT(L) 

IF (TID-TIOH) 10,9,11 
J = J + l 

CFS( JI=DCFS(I ,L) 

T IO=TID>Or (K) 

GO TO 11 
J = J+1 

^CFS(J)=OCFS(I-lfLHI(T IO--T lOH+OT (L) I /DT( L) l^(OCFS( I , L )-nCFS ( I - I, L > 

TIC=na+DT(K> 

CONTINUE 

IEND(L)=J 

DT(L)^DT{K) 

DO 12 1=2, J 
UCFS(I,L»=CFS( n 

IF ( IEND( I DU-IENO( 1021) 14,14,15 
M= IEND( IDl ) 

GO TO 16 
M= IEND( ID2 J 
IF (M'-MAX) 17,17,18 
DETERMINE TIME SCALE 
MRTQ=MAX/M 
XMRTO=MRTO 
GO TO 19 
M = MAX 

YSCL=QM AX/50. 

PLOT HYDROGRAPHS 
on 20 1=1, MAX 
CFS( I )=DASH 

WRITE (6,41) QMAX,(CFS( 1) ,I = l,MAX),OOT 

Q1=QMAX 

Jl=lO 

00 37 J=I,50 

IF (J-Jl) 23,21,23 

DO 22 1=1, MAX 

CFS( I)=DASH 

GO TO 25 

00 24 1=1, MAX 

CFSI I) =BLANK 

Q2=Q1-YSCL 

K=1 

DO 28 1=2, M 
K=K+MRTO 

IF (OCFSI I ,IDI)"QU 26,27,28 
IF (OCFSI I , IDl )*-Q2) 28,28,27 
CFS(K)=ZERO 
CONTINUE 

WRITE (6,44) DOT, (CFSI I) ,1=1, MAX) , DOT 

IF (102) 34,34,29 

K=l 

DO 33 1=2, M 


T 26 
I 27 
I 28 
r 29 
I 30 
I 31 
I 32 
I 33 
I 34 
t 35 
I 36 
I 37 
I 38 
I 39 
I 40 
1 41 
I 42 
I 43 
I 44 
1 45 
1 46 
I 47 
1 48 
I 49 
I 50 
I 51 
I 52 
I 53 
I 54 
I 55 
I 56 
I 57 
I 5B 
I 59 
I 60 
I 61 
I 62 
I 63 
1 64 
I 65 
I 66 
I 67 
I 68 
1 69 
I 70 
I 71 
I 72 
I 73 
I 74 
I 75 
I 76 
I 77 
I 78 
I 79 
I BO 
I 81 
I 82 
I 83 
I 04 
I 85 


63 



461 


462 


463 

30 

464 

31 

465 


466 

32 

46 7 

33 

46 8 


46 9 

34 

47 0 

35 

471 


472 

36 

47 3 

37 

474 


475 

C 

476 


477 

38 

478 


479 


4 80 


481 

39 

482 

40 


C 

483 

41 

484 

42 

485 

43 

48 6 

44 

487 

45 

488 

46 

409 

47 

490 


491 

C 

492 



493 


494 


495 


496 


497 

C 

498 


499 

1 

500 


50 1 


502 


503 

2 

504 


505 


506 


507 

3 

508 


509 

4 

510 


511 


512 



K = K+MRTO 

IF (OCFS( I #ID2)-Q1) 30t31,32 

IF (OCFS( I tID2»-Q2» 32,32,31 

CFS(K)=PLUS 

GO TO 33 

CFS(Ki==BLANK 

CONTINUE 

WRITE (6,42) (CFS( I) ,I=1,MAX) 

IF (J-Jl) 36,35,36 

JUJl+lO 

WRITE (6,43) Q2 

Ql=Q2 

CONTINUE 

CFS( UPTIME 

OTT=OT( 101)^10. /XMRTO 

PUT TIME ARRAY IN CFS AND WRITE TIME SCALE 

00 38 1=^2, 12 

CFS( n-CFSC I-U+OTT 

WRITE (6,45) (CFSd ), I^'lfU) 

WRITE (6,46) 

IF (NPU) 40,40,39 
WRITE (7,47) 101,102 
RETURN 

FaRMAT(lX,F640,« CFS.SU9A1) 

FORMAT ( lH+,llX,liaAl) 

FORMAT { IH+,F6*0) 

FORMAT (11X,120AI) 

FORMAT (6X,12F10#3) 

F0RMAT(49X,*TIME HOURSV//) 

FORMAT! *PLOT HYDST2lfMD 1 = ' , II ,T29, M D II^SIl) 
END 


I 86 
I 87 
I 88 
I 89 
I 90 
I 91 
I 92 
I 93 
I 94 
I 95 
I 96 
I 97 
I 98 
I 99 
I 100 
I 101 
I 102 
I 103 
I 104 
I 105 
1 106 
I 107 
I 100 
I 109 
I 110 
I 111 
I 112 
I 113 
I 114 
I U5 
I 116 
I U7- 


SUBROUTINE ADHYD J I 

THIS SUBROUTINE ADDS TWO HYDROGRAPHS# J 2 

COMMON CFS (300) ,OCFS( 300, 6) , lEND ( 6) ,DATA( 310) , DA(6 ) , DP( 20 ) , NPU ,NHD J 3 

l,SCFS(20),C(20) ,A(20,6),Q(20,6),RAIN(200) , DEEP ( 20 , 6 ) , NER , MAXNO , NCO J 4 

2MM,ICC,NCO0E,DIST(6),SEGN(6),CTBLE(50,ll) , ITBLE ( 50 ,2 ) ,2ALFA( 20 ) ,0T J 5 

3(6), TIME, PEAK(6),R0INtISG(6) J 6 

ID=DATA(1I J 7 

NHD=DATA(2) U 8 

ID1*0ATA(3) J 9 

ID2«0ATA(4) J 10 

PEAK(10)-1. J U 

MAKE TIME INCREMENTS EQUAL IF NOT EQUAL# USE SMALLER INCREMENT J 12 

IF {0T( IDl )-DT( ID2) ) 1,3,2 J 13 

OT(IO)^Or( lOl ) J 14 

L^lDl J 15 

K==ID2 J 16 

GO TO 6 J 17 

OT(ID)=OT( 102) J 10 

L=:ID2 J 19 

K=IOl J 20 

GO TO 6 J 21 

0T( ID)«DT( IDl) J 22 

IF (IEN0(ID1)-IEND(I02)) 4,4,5 J 23 

M3-IEND(1DI) J 24 

KUID2 sJ 25 

IEND(I0)=IEND(I02) J 26 

GO TO 18 J 27 


64 



*^13 

5 

514 


515 


516 

C 

5L7 

6 

518 


5L9 


520 


521 


522 

7 

523 


524 


52 3 


52 6 

8 

527 


52 8 


52 9 

9 

530 

10 

531 

11 

532 

C 

533 


534 


535 


536 


537 

12 

538 

13 

539 


540 


541 


542 

14 

543 


544 


545 


546 

15 

54 7 

16 

54 8 


.549 

17 

550 

18 

551 


552 

C 

553 


554 


555 


556 

19 

55 7 

20 

55 8 


559 

21 

560 

22 

561 

23 

562 


563 


564 

24 

565 

25 

566 


567 

26 

568 

27 


M3=I END( ID2) 

K1=ID1 J 28 

IENO( IOI=IENO< lOU J 29 

GO TO 18 0 30 

DETERMINE DURATIONS OF FLOW J 31 

X I END1= lENOI IOl)-l J 32 

X1EN02=IEND(ID21-1 J 33 

DUR1 = X lENDl’l'OTI IDU J 3X> 

DUR2=XIEND2>l'DT( ID2) J 35 

IF (DUR1-DUR2) 7,8,8 J 36 

lENDl I0|^0UR2/DT< lOlFl. J 37 

H3=DUR1/0T( lOI+l. J 38 

K 1=1 02 J 39 

GO TO 9 J ^0 

IEN0(IDI=0UR1/DT(IDU1. J 

M3=DUR2/DT(ID» + 1. J '•a 

K1=ID1 J 43 

IF ( lENOI ID)-3001 11,11,10 

IEND(ID) = 300 '> 45 

M2 = ieNr)(K) J 46 

J B 1 J 47 

interpolate ONE HYDROGRAPH IF NECESSARY 

TIDH=0. J 49 

TI0 = DT( ID) J 50 

00 15 1=2, M2 J 51 

TIDH=T IDHFDTI K) J 52 

IF (TIDH-TIO) 15,13,14 ^ 53 

J=J+1 J 54 

OATA ( J)=0CFS(I,K» 55 

TID=TIDfOT( ID) -J 56 

IF (J-3001 15,16, 16 57 

J=J+1 J 58 


iDATA <J>=nCFSn-l,K) + ((TID-TIDHFDT(K))/OT(K)»*(OCFS(I,K)-OCFS(I-l, J 60 

TID=TIO+DT( 10) 

IF (J-300) 12,16,16 

continue 

IEND(K)=J 
DO 17 1=2, J 
OCFSI l,K)=DATA( I) 

M=IEND( ID) 

DAIID) = DA( ID1)F0A<ID2) 

R0=0. 

ADD HYDROGRAPHS 
DO 20 1=1, H3 

OCFSI I , ID)=OCFS( I , IDl )+OCFS( I , ID2) 

IF (OCFSI I , ID)-PEAK(IO) 1 20,20,19 
PEAK(ID)=OCFS(I,IDI 
RO=RO+OCFSI 1,10) 

IF (PEAK! IC)-PEAKIK1) I 21,22,22 
PEAKI ID1=PEAK(K1I 
IF IM-M3) 25,25,23 
M3=M3tl 
DO 24 I=M3,M 
OCFSI I , ID)=OCFSI I ,K1) 

RO=RO+OCFS( I , ID) 

ROIN=I RO*OT( I D) ) / lOAI IDl *64 5,333 ) 

IF (NPU) 27,27,26 
WRITE (7,28) ID,NHO r 101 , 102 
RETURN 


J 61 
J 62 
J 63 
J 64 
J 65 
J 66 
J 67 
J 68 
J 69 
J 70 
J 71 
J 72 
J 73 
J 74 
J 75 
J 76 
J 77 
J 78 
J 79 
J 80 
J 81 
J 82 
J 83 
J 84 
J 85 
J 86 
J 87 


65 




c 


J 

88 

569 

28 

FORMATl »A00 HYO» rT2lf MD==* f 11 1729, • HYD N0== • * I 3 , T45 f » I D 

J 

89 



IT60,»ID iI = Ml) 

J 

90 

570 


END 

J 

91- 

571 


SUBROUTINE SRC 

K 

l 


C 

THIS SUBROUTINE STORES AN ELEVATION - END AREA - FLOW TABLE. 

K 

2 

572 


COMMON CFS (300) fOCFSI 300,6) , I END (6) , DATA ( 3 10 ) , OA ( 6 ) ,DP( 20) ,NPU,NHD 

K 

3 



1, SCFS( 20 1 ,C (20) f At 20 ,6) ,Q(20f6 ) ♦ RAIN (200 ) f DEEP (20, 6) , NER , MAXNO , NCO 

K 

4 



2MM ,ICC,NCQDE,DIST(6) ,SEGN(6) ,CTBLE (50 ,11) , I TBL E ( 50 , 2 ) , ZAL F A ( 20 ) , OT 

K 

5 



3(6), TIH6fPEAK(6) ,R0IN,ISG(6) 

K 

6 

573 


10=0ATA(1) 

K 

7 

574 


VS=DATA( 2) 

K 

B 


C 

VALLEY SECTION NUMBER 

K 

9 


c 

REMAINING DATA ARE ELEVATION, AREA, AND FLOW FOR EACH POINT OF 

K 

10 


c 

THE rating CURVE 

K 

11 

575 


EHIN^DATAO) 

K 

12 

576 


J==3 

K 

13 

577 


DO 1 1^1,20 

K 

14 

578 


ELEV=OATA( J) 

K 

15 

579 


DEEP( I ,IO)=OATA( J)-EMIN 

K 

16 

580 


A( I ,IO}=DATA( J + l) 

K 

17 

581 


Q{ I,ID) = DATA( J+2) 

K 

18 

582 


J^J+3 

K 

19 

583 

1 

CONTINUE 

K 

20 

584 


RETURN 

K 

21 

585 


END 

K 

2 2- 

586 


SUBROUTINE CMPKC 

L 

1 


c 

THIS SUBROUTINE COMPUTES THE DISCHARGE END-AREA ELEVATION 

L 

2 


c 

RELATIONSHIP FOR A VALLEY SECTION. 

1 

3 

587 


COMMON CFS (300) ,OCFS( 300,6) , IEND(6) , DATA ( 3 10 ) , DA ( 6 1 , 0 P ( 20 ) ,NPU,NHD 

L 

4 



l,SCFS{ 20) ,C(20) ,A(20,6) ,Q (20, 6) , RA1N( 200 ) , DEEP ( 20 , 6 ) ,NER , MAXNO, NCO 

L 

5 



2MM,ICC,NC00E,DIST(6),SEGN{6),CTBLE(50,ll) , I TBL E ( 50 , 2 ) , Z ALF A ( 20 ) , DT 

L 

6 



3( 6) ,TIME,PeAK(6) ,R0 IN, ISG(6) 

L 

7 

588 


ID==OATA( 1) 

L 

a 


c 

STORAGE LOCATION NUMBER. (1-6) 

L 

9 

589 


VS=0ATA(2) 

L 

10 


c 

VALLEY SECTION IDENTIFICATION NUMBER. 

L 

LI 

590 


NSEG^OATAt 3) 

L 

12 


c 

NUMBER OF SEGMENTS IN THE VALLEY SECTION. 

L 

13 

591 


ELC==0ATA(4) 

L 

14 

592 


EMAX^DATA{ 5) 

L 

15 


c 

MAXIMUM elevation FOR COMPUTATIONS. 

L 

16 

593 


SL0PE1=DATA(6) 

L 

17 


c 

CHANNEL SLOPE. 

L 

18 

594 


SLOPE2=DATA(7) 

L 

19 


c 

FLOOD PLAIN SLOPE. 

L 

20 

59 5 


DIF=(EMAX-EL0)/19. 

L 

21 

596 


C{ 1)=EL0 

L 

22 

597 


DO 1 1=^2, 20 

L 

23 

598 

1 

C( n=C( I-l )+DIF 

L 

24 


c 

SET AREA AND DISCHARGE ARRAYS = 0. 

L 

25 

599 


DO 2 1=1,20 

L 

26 

600 


A( 1,101=0. 

L 

2 7 

601 

2 

Q{ I,ID)=0. 

L 

28 

602 


J = 8 

L 

2 9 

60 3 


WRITE (6,24) VS 

L 

30 


c 

READ N VALUES AND SEGMENT BORDER POINTS. 

L 

3L 

604 


DO 3 r=l,NSEG 

L 

32 


66 



605 


606 


607 

3 


C 

608 


609 


610 

6 

611 


612 

5 

613 

6 


C 

614 


615 


616 


617 


618 

7 

619 


620 


621 

8 

622 

9 


C 

623 


626 


625 


626 


62 7 


62 8 

10 

629 


630 

11 

631 

12 

632 

13 

633 


636 

16 

63 5 


636 


637 

15 

636 


639 

16 

660 


661 


662 

17 

66 3 

18 

646 


665 


666 

19 

66 7 


66 8 


669 


650 


651 


652 


653 


656 

20 

655 



c 

c 

656 

657 

658 21 

659 


SEGN( I ) =OATA( J J 
D !ST( n:=^OATA( J+1 ) 
J = j4-2 


REJAIN-[NG DATA ITEMS ARE DISTANCES AND ELEVATIONS. 

DO 6 I=1,NSEG 
J = J + 2 


IF (DATA! J)-OIST(I| I A, 5, 5 

ISGI l)=J+l 

CONTINUE 

COMPUTE DISCHARGES AND END 
DO 22 K=1,NSEG 
J* JJ J 

JJ 1 
IF (SEGNIKI) 7,7,0 
SLOPe=SLOPEl 
SEGNIK ) = -SEGN ( K ) 

GO TO 9 


areas FDR BACH SEGMENT. 


SL CPE=SLf)PF:2 

SLPN = l.486’>SL0PE’l‘’*.5 

COMPUTE AREA AND DISCHARGE FOR SEGMENT. 

UU <:J. I-2t20 

AA=0. 


P = 0. 


J=JJJ-1 

DEP2=0. 

J = J + 2 

IF (J-ISG(K)) 12,12,11 
IF (AA-,OOn 21,21,20 
IF (DATA! J)-C< I ) I 13, 10,10 
OEPl=Cn )-DATA( J) 

IF (J"JJJ1| 16,16,1a 
XL=DATA( J-1 )-DATA( J~3) 

DEP3=ABS(OATA(J-2)-DATA(J>) 

XL=XL’«DEP1 /DEP3 
AA=AA+XC*( OEP 1+06P2 )/2. 

P=P+SQRT ( ( OEP 1-DEP2)*>('2*XL'«^*2I 
DEP2=DEP1 
J = J+2 

IF (J-ISG(K)I 17,17,20 
IF IDATAIJ |-C( I M IS, 10,19 
DEP1=CI I )-DATA(J) 

XL=DATAIJ-1)-0ATAIJ-3| 

GO TO 15 
0EP1=0. 

XL=DATA( J-1)-DATA( J-3> 

DEP3=ABS(OATA(J-2 »-OATA(J I ) 

XL=XL*DEP2/DEP3 

AA=AA+XL*(DEP1+DEP2)/2. 

P = P+SQRT< ( DEPl-DEP2)’t'*2+XL**2l 
0EP2=0. 

GO TO 10 
R=AA/P 

SGN=SEGNI K 0025*R 

^ n™!" totals fo« 

m' 1 : IS r.Jl !; 1 s 1 

CaNTINU(= 


JJ J=J-3 


L 33 
L 36 
L 35 
L 36 
L 37 
L 38 
L 39 
L 60 
L 61 
L 62 
L 63 
L 66 
L 65 
L 66 
L 67 
L 68 
L 69 
L 50 
L 51 
L 52 
L 53 
I 56 
L 55 
L 56 
L 57 
L 58 
L 59 
L 60 
L 61 
L 62 
L 63 
L 66 
L 65 
L 66 
L 67 
L 68 
L 69 
L 70 
L 71 
L 72 
L 73 
L 76 
L 75 
t 76 
L 77 
L 78 
L 79 
L 80 
L 81 
L 82 
L 83 
L 84 
L 85 
L 86 
L 87 
L 88 
L 89 
L 90 
L 91 
L 92 


67 



660 

22 

CONTINUE 

L 

93 

661 


DO 23 I=l»20 

1 

94 

662 


oeEP(i,io)-c(n-ELO 

L 

95 

663 


WRITE (6,25) C( 1 » tA( I , ID) fQ< If ID) 

L 

96 

664 

23 

CONTINUE 

L 

97 

665 


RETURN 

L 

98 


c 


L 

99 

666 

24 

FORMAT ( LHO fT42t CRATING CURVE VALLEY SECT ION •, F 5 . 1/ T46 ,« WATER *» T56 , 

L 

LOO 



PFLQW* ,T66,*FL0WVT45f» SURFACE* ,T56,»AREA' , T66 , • R ATE • /T46 , • EL E V * tT 

L 

LOl 



256,* SQ FT* »T66, *CFS* ) 

L 

102 

667 

25 

FORMAT (40XtF l0*2f2F10.1) 

L 

103 

668 


END 

L 

104- 

669 


SUBROUTINE STT 

M 

1 


C 

THIS SUBROUTINE STORES A DEPTH - FLOW - TRAVEL TIME TABLE. 

M 

2 

670 


COMMON CFS (300) fOCFSt 300,6) , I eND|6 ) ,0ATA(310) , DA ( 6 ) , DP ( 20 ) , NPU fNHD 

M 

3 



1, SCFS{ 20) tC(20) ,A(20,6l , Q ( 20 1 6 ) , R A I N ( 200 ) , DEEP ( 20 , 6 ), NER , MA XNO , NCO 

M 

4 



2MM,ICC,NCODE,OIST{6) , SEGN <6 ) ,CTBLE ( 50 , 1 1) , I TBL E( 50 , 2 1 , Z ALF A ( 20 ) , DT 

M 

5 



3( 6),TIHE,PEAK(6) fROIN, ISG(6) 

M 

6 

671 


IO=DATA( 1) 

M 

7 

672 


REACH= DATA (2) 

M 

8 

673 


XL=DATA(3) 

M 

9 

6 74 


SLQPE=0ATA(4) 

M 

10 

675 


DIST( ID)=SLQPe*XL 

M 

n 

676 


J^5 

M 

12 

677 


00 1 I==l,19 

M 

13 

678 


OP ( I) = DATA( J) 

M 

14 

679 


SCFSd )=DATA( J+1) 

M 

15 

680 


C( II^OATAI J+2) 

M 

16 

681 

1 

J = J+3 

M 

17 

682 


RETURN 

M 

18 

683 


END 

M 

19- 

6 84 


SUBROUTINE CMPTT 

N 

1 


C 

THIS SUBROUTINE COMPUTES THE TRAVEL TIME AT GIVEN 

N 

2 


c 

DISCHARGE RATES 

N 

3 

635 


COMMON CFS (300) fOCFSI 300,6) , IEND{6) , DATA ( 3 10 ) , DA ( 6 ) , OP ( 20 ) , NPU ,NHO 

N 

4 



1,SCFS(20| ,C(20) ,A(20,6),Q(20,6) , RA1N( 200 ) , DEEP ( 20 , 6 ) , NER, MAXNO,NCO 

N 

5 



2MM,ICC,NC0DE,DIST(6) ,SEGN(6) ,CTBLE (50, lU , ITBLE( 50,2) ,ZALFA( 20 ) , DT 

N 

6 



3( 6),TIME,PEAK(6) ,ROIN, ISG(6) 

N 

7 

686 


IO=DATA( 1) 

N 

8 

687 


ReACH=DATA(2) 

N 

9 

688 


N0VS=DATA(3) 

N 

10 

689 


XL^0ATA(4) 

N 

11 

690 


SL0PE=DATA(5) 

N 

12 

691 


DIST(IO)=SLQPE*XL 

N 

13 

692 


XLD36=XL/3600. 

N 

14 


c 

ZERO ARRAYS 

N 

15 

693 


DO 1 J=l,20 

N 

16 

694 


DATA (J)=0. 

N 

17 

695 

1 

CFS( J)=0. 

N 

18 

696 


ID1=1 

N 

19 


c 

FIND RATING CURVE WITH SMALLEST MAXIMUM FLOW RATE 

N 

20 

697 

2 

QMIN=Q(20, lOl) 

N 

21 

698 


MIN^IOl 

N 

22 

699 


GO TO 4 

N 

23 

700 

3 

IDl=lDlfl 

N 

24 

701 


IF (QMIN-Q(20,ID1)) 4,4,2 

N 

25 

702 

4 

IF (lOl-NOVS) 3,5,5 

N 

26 

703 

5 

I^l 

N 

27 



7Q5 

706 6 
C 

707 
700 
709 

71 0 

111 1 
7L2 L7 
7L3 
7L^ 

715 8 

716 

717 

718 9 

719 

72 0 

C 

721 
72 2 
72 3 
724 
72 5 

72 6 

727 10 

C 

726 

729 11 

73 0 

73 1 12 

C 

732 13 


733 14 

734 15 


73 5 16 

736 

737 

C 

C 

c 

738 


739 
74 0 

741 

742 

743 

744 

C 

74 5 

746 I 


^ 'X'"™ CURVE 

SCFS( I (=g( J,MIN) 

I =t+l 

CUMPUT ENIi AREA AND DEPTH 

00 9 idi=i,novs 

DO 9 J=l,19 
DO 7 I =2,20 

TF (0( [ , IDl )-SCFS( J ) » 7,17,8 
CONTINUE 

DATA ( J)=A( I, I01)+DATAt J) 

CFSI J»=()EEP( I ,ID1 ItCFSCJ) 

GO TO 9 

XY = ( SCFSI J )-gn-l,IDin/( 0( I , IDll-Q{ I-I mill 

XNCVS=NOVS 
WRITE (6,13) REACH 
COMPUTE TRAVEL TIME 
DO 10 1=1,19 
AVAREA=OATAI I )/XNOVS 
OP(n = CFS( I J/XNOVS 
S= AVAREA*XL036 
Cl I l=S/SCFS( I J 

WRITE (6,1A) 0P( I ) ,SCFS( I ) ,C( I ) 

CONTINUE 

PUNCH CODE 

IF (NPU) 12,12,11 

WRITE (7,15) ID, REACH, XL, SLOPE 

fO'’<n.SCFS(n,C(n,I=l,19) 

2',T56, •CFS'.T66,'HRS*) • T56, ‘RATE* , T66 ,' T IME •/ TA6 ,' FEE 

FORMAT ('tOX,F10.2,F10.0,F10.4) 

l'LENGTl=s??!oI”^??vJa!!sLOPE=SF8’lIlJ?/Fn'^*'^^ NQ= • ,f 5. 1 , r44 , 
END 


suBRourmE route 

I 0=DATA( 1 > 

NH0=DATA(2 ) 

IDH=DATA(3) 

0TnD)=DATA(4) 

DA(ID)=f)A< IDH) 

M= I END! IDH) 

IF ID AND IDH ARE EQUAL, ADD I TO IDH 

IF no- IDH) 3,1,3 

IOh=IDH+l 


N 28 
N 29 
N 30 
N 31 
N 32 
N 33 
N 3A 
N 35 
N 36 
N 37 
N 38 
N 39 
N AO 
N Al 
N A2 
N A3 
N AA 
N 45 
N 46 
N 47 
N 48 
N 49 
N 50 
N 51 
N 52 
N 53 
N 54 
N 55 
N 56 
N 57 
N 58 
N 59 
N 60 
N 61 
N 62 
N 63 
N 64 
N 65 
N 66 
N 67 
N 68 
N 69- 

0 1 

0 2 
0 3 

□ 4 
D 5 
0 6 
D 7 
0 8 
O 9 

□ 10 

0 11 
0 12 
□ 13 

0 14 

0 15 

0 16 
Q 17 


69 



74 7 


74 8 

2 

74 9 


750 


751 

3 

752 


753 


754 


755 


756 


757 


758 


759 


760 

C 


C 

76 1 


762 

4 

763 


76 4 


76 5 


766 


76 7 

5 

76 8 


76 9 


770 


771 

6 

77 2 


773 


774 

7 

775 


776 

8 

777 


77 8 


77 9 


780 

9 

781 

10 

782 


78 3 


784 


7R5 

il 

78 6 


787 


78 8 


789 

12 

790 

13 

79 1 


792 


793 


794 

14 


C 

795 

15 

796 


797 

16 


C 

796 

49 

799 
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DO 2 1=1, M 

OC FS< I ,n)H)=OCFS( I, tDH-1) 

DTI IDHI=DT(IOH-l) 

PEAK! I DH( = PEAK( lDH-1) 

NePRT=0 
PE AKI I D) =1 . 

RO=0. 

N=19 

OCFSd ,IO)=0. 

S = 0. 

Ti=cm 

j=i 

GUES=l. 

CFSIl) =0. 

IF ROUTING INTERVAL IS NOT EQUAL TO TIME INCREMENT OF INFLOW 
HYDROGRAPH, INTERPOLATE 
IF (DT{ID»-OT( IDHI) 8fl5,<i 
T I D= DT (ID) 

T I 0H=0 . 

DO 7 1=2, M 
TIOH=T 1DH+0T( lOH) 

IF (TID-TIOH) 6,5,7 
J = J+ 1 

CFSI J ) =OCFSI I I lOH) 

TIC=TIO+DT(ID) 

GO TO 7 
J = J+1 

CFSI J)=OCFS( I - I , IDH» + ( ( TID-T lOH+OT (IDHH /DTII DH M * ( OCf S ( I , I DH)-OC I 
1S( I-l, IDHI ) 

TIC=TID+DTIID) 

CONTINUE 
GO TO 13 
T I 0H=0 . 

T I D=OT IIDI 
DO 12 1=2, M 
T I DH=T IDH+DT I IDHI 
IF (TIDH-TIDI 12,10,11 
J = J<-1 

CF S( J) =OCf SI I , lOH) 

T I D=TI D+DTI ID J 
IF IJ-3001 12,13,13 
J= J+1 

CFSI Jl =OCFSI I-l, IOH) + I (TIO-TIDH+OTIIDHII/OTI IDHI l*IOCFSI I , IOH)-OC I 
ISI I-l, IDHI I 
TIC=TID + 0Tn01 
IF IJ-300) 9,13,13 
CONTINUE 
lENDl IOH>=J 
DT(1DH)=UTI10I 
H = J 

00 19 1=2, M 

OCFSI I , IDH)=CFSII I 

IF INFLOW IS ZERO, SO IS OUTFLOW 

DO 16 L=2,M 

IF tOCFSIL, IDHI I 16,16,99 
OCFSIL ,ID1=0. 

ROUTE 

DATA IL-1I=0. 

DO 92 I=L,300 
IF (I-M) 18,18,17 



801 

17 

802 

18 

803 


804 

C 

805 


806 

19 

807 


806 


809 

20 

810 


fill 

21 

812 


813 

22 

814 


815 


816 

23 

817 


818 


819 

24 

820 


021 


822 

25 

823 

C 

824 


82 5 

26 

82 6 


82 7 


82 8 

27 

82 9 


830 

28 

831 


832 

29 

833 


834 


835 

30 

836 


837 

C 

838 

31 

839 


840 

32 

841 


842 

33 

843 


844 

C 

845 


84 6 


04 7 


84 8 


849 

C 

85 0 


851 

34 

852 

35 

853 


854 


855 



oc FS( r ♦ I OH ) -DCFS ( I*- I # I DH ) ^8^ , 9 

AVIN-(DCFSn, IDH)fOCFS(I-l,iDHn/2 

S I A=fS + AVIN ^ 

J^l 

DETERMINE DEPTH AND TRAVEL TIMF np r. 

IF lOCFSd.rOH.-SCFSn?) I9l2;f20 

I • I OH I /SCFS ( 1 ) ) ♦OP ( 1 ) 

GO TO 25 
DO 21 J=2,N 

IF < OCF S ( I I 1 DHI-SCFS( J 1 ) 24t23.21 
CONTINUE 

IF (NERRT) 22,22,36 
WRITE (6,A6) 

NE RRT=1 
GO TO 36 
D I 2- DP ( J » 

TI2=C( Jl 
GO TO 25 

RATIO=(OCFS(I,IOH»-SCFS(J-l))/r<:rPC( .. cr... 
DI2=DP(J-l)+RAT10*t0P(j -DP j^[n 

TI2=C<J-1)^RATI0*{C(J)-C(J-- H 

00 35 IT=l,IO ' 

002*=IGUES/SCFS< 1) )*DP{ U 
T0 2 = C(n 
GO TO 31 
DO 2fl J = 2,N 

IF <GUBS-SCFS(J)> 30,29,28 
CONTINUE 
J = N 


D02=DP( J) 
T02=C( Jl 
GO TO 31 


R A Tla= I GUR S-SCFS( J-l I I / ( SCFSIJI-SCF 51 1-11 i 
D02-DP ( J-l l+R AT lO+l DP{ J)-OP( j_i J ) 

T0 2=C( J-l) +RATIO*(C(J)-C(J-ni 
FIND WATER SURFACE SLOPE 
DUD=DISTI ID)/(OIST( 101+012-002) 

IF (ODD-. 01) 32,32,33 
GUeS=OCFS( I-l , lOH) 


GO TO 35 
T2=.5*<TI2+T02| 
T2=T2>KSQRT(DDD) 
T=TI+T2 


CON-PUTE ROUTING COEFFICIENT 
C0EF = (2.«f)T( IO))/(r+DT( ID)) 
02=C0EFi‘Sl A 
rRYl=GUES 

R ATI 0=02/ I CUE S+, IE- 20 I 
OIFF=AHS( 1. -RATIO) 

TEST FUR CONVERGENCE 
IF (OIFF-.OOI) 37,37,34 
GUES=n 2 
CONTINUE 

OCFS( I,II))=:OATAn-l)*SIA 

DATA ( I )=|JATA( I-U 
WRITE (6,471 1,0CFS(I,ID) 
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0 79 

0 80 
0 81 
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0 84 

0 05 

□ 86 
0 87 

0 88 
0 89 

0 90 

0 91 

0 92 

0 93 

0 94 

0 95 

0 96 

0 97 

0 90 

0 99 

0 ion 
0 ini 
n 102 

0 103 
0 104 
0 105 
0 106 
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0 IPB 
0 109 
0 110 
0 111 
0 112 
0 113 
0 114 
0 115 
0 116 
0 117 
0 lie 
0 119 
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n 129 
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0 135 
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856 


857 

36 

858 


859 


86 0 

37 

861 

C 

862 

38 

863 


864 


865 


866 

39 

86 7 

40 

868 

41 

869 

42 

870 


871 

43 

872 

C 

873 


874 

44 

875 

45 


C 

8 76 

46 

877 

47 

878 

48 

87 9 


88 0 

C 


C 

881 



882 


883 


864 


805 


886 


887 


888 


889 


890 


891 

C 

892 


893 


894 

C 

895 

1 

896 


897 


898 


899 

2 

900 


901 


902 

3 


GO TO 38 

QCFS{ l,IO)=DATA{ I-1»*SIA 
DATA ( I)=OATA(I-l) 

GO TO 38 
0CFSII,I0)=02 
DATA in=COEF 
COKMUTE NEW STORAGE 
S=SiA-OCfS(I,IDI 
T1=T2 

RO-RO+OCFS( 1 1 ID» 

IF (OCFS( 1 ,ID)-OCFS( 1-lfIOn 39fA0,A0 
IF (OCFSt 1,ID>-1.) A3,43,A2 
IF (OCFSt I ,IO)*PBAK(ID) ) 42,42,41 
PEAK(IOI=OCFStI,ID) 

CONTINUE 
1 = 300 

I £ NO t I D I = I 

ROIN=(RO*OT( 101 »/(DA{ 101*645.333) 

PUNCH CODE 

IF (NPUl 45,45,44 

WRITE (7,481 ID,NHD,IOH,DT( ID) 

RETURN 

FORMAT! IHO, 'TRAVEL TIME TABLE EXCEEDED') 

FORMATITIO, 'PROBLEM FAILED TO CONVERGE AFTERIO ITERATIONS. CONVERG 
lENCE WAS FORCED. '/TaO, 'OUTFLOW NUMBER = ',14, 'RATE =',F10,2> 
FORMAT! 'ROUTE', T21, '10 = ', II, T29, 'HYD N0=*, 13, T45, 'INFLOW ID=',I 
11 ,T65,'DT=',F8.6,'HRS'I 
END 


0 138 
0 139 
0 140 
0 141 
n 142 
0 143 
0 144 
0 145 
0 146 
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0 149 
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0 160 
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0 163 
0 164 
0 165 
0 166- 


SUBROUriNE RESVQ P 1 

THIS SUBROUTINE ROUTES A HYDROGRAPH THROUGH A RESERVOIR WITH THE P 2 

STORAGE-INDICATION METHOD. P 3 

COMMON CFS( 300) ,OCFS(300,6) , 1END(6),0ATA(310) ,DA(6) ,DPI 20) ,NPU,NHD P 4 

1, SCFS( 20) ,C{ 20) ,A(20,6) ,Q( 20 ,6 ) , RA I N( 200 ) , DE EP ( 20 , 6 ) , NER , MA XNO, NCO P 5 

2MM,ICC,NCOOE,DIST(6),SEGN(6) ,CTBLE ( 50 , 1 1 ) , I TBLE ( 50 , 2 ) , ZALF A ( 20 ) , DT P 6 

3(6), TIME, PEAK(6),R0IN,1SG(6) P 7 

ID=DATA(l) P S 

NH0=DATA(2) P 9 

1DH=DATA(3) P 10 

NERES=0 P 11 

0T(1D)=DT| IDH) P 12 

R0=0. P 13 

DA ( ID)=0A( IDH) P 14 

PEAKI1D)=1. P 15 

J=1 P 16 

I»4 P 17 

REMAINING DATA ARE FLOW AND STORAGE VALUES P 18 

SCFS(J)=DATA( I) P 19 

STCRE1=0ATA(I+1)*12.1 P 20 

ST0RE=ST0RE1 P 21 

COMPUTE STORAGE COEFFICIENT ARRAY C P 22 

C( Jl=(SCFS(J)/2.)+(ST0REZDT(ID)) P 23 

I = H-2 P 24 

J=J+1 P 25 

IF (J-20) 2,2,3 P 26 

SCFS(J)=DATA( I) P 27 

ST0RE=DATA(1+1)*12.1 P 28 

IF (SCFS( J)-.001) 3,3,1 P 29 

N=J-1 P 30 
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903 

904 

C 

905 

906 

907 4 

908 5 

909 

C 

910 

911 

912 6 

913 

914 7 

915 

916 8 
C 

917 
910 

919 9 

920 

921 10 
C 

922 1 1 
92 3 

92 4 

925 12 

926 13 

927 14 

928 15 
92 9 

930 16 

931 

C 

932 

933 17 

934 

935 

936 10 
C 

937 19 

938 20 


939 21 

940 

941 

C 

C 

942 


943 

944 

945 

946 

947 


OCf-S( 1 , U)| =0. 

s=srniui/r)T(ic) 

Kourc 

on 15 1=2,150 

11^ ( I-IEND( ION) ) 5,5,4 

t)CFS{ I , IQH)=0.0 

AVIN=(GCFS{ I, IDHlfDCFSd 

S 1 A=S+AVIN 

DETCKMINH PROPCK C 

DO 6 J=1,N 

IF (.SIA-CfJ)) 10,9,6 

COMTINUr 


“It IDHI 1/2. 


IF (MFIM'S) 7,7,8 
WfUTK (6,19) 
NF9FS=1 


R F SC = SC(-S ( N ) /C<N ) 
CfJHPUT OUTFLOW 
Of.FSI I , ID) =RnSC*SIA 
GO TCI 11 


OCFSI I , ID)=SCFS( Jl 
GO TO 11 


nCFSI I , I0)=SCFSIJ- 
11 ) ) 


1 >+( (SIA-C(J- 


n)/(C(J)-C(J-l)))*(SCFS(J)-SCFS(J- 


OFTRRMINF NEW STORAGE 
S = S] A-OC)=S( I, ID) 
f<0=IU)*'0CFS( 1,11)) 


F KJCFSl I ,I1))-0C)-SII-1,|D)) 12,13,13 
It- (OCFS(I,I())-l.) 16,16,15 
IF (UCFSl I , ID)-PEAK(1DI I 15,15,1/, 
l'EAK(II)) = OCFS(I,II)| 

CONT INUE 


I = 150 


lFNI)(ID)cI 

Rl)IN = Rn-l<l)TnO)/(DA( II))«6/,5.333) 
PUNCH CODE 
IF (NPU) 18,10,17 
WRITE (7,20) ID,N)I(),II)H 
I I «2»N + 3 

WRITE (7,2 1) ll)ATA( I ) , Ia5,i I I 
RETURN 


FORMAT ( IHO , 33IISTr)RAGE-OI SCHARGE 
FORMAT( •ROUTE RE SER VOl R * , T?1 i 
ILUW 11)=', II 
2A(;E(AC FT)*) 

FORMAT (T2l,F10.l,F13.l) 

END 


table exceeded,) 

ID=' I Il,T29, 'HYD N0» • , 1 3 , TA2 , * INF 
/T21,»0UTFL0W(CFSI',T37,'ST0R 


P 31 
P 32 
P 33 
P 3A 
P 35 
P 36 
P 37 
P 38 
P 39 
P AO 
P A1 
P 62 
P 43 
P 44 
P 45 
P 46 
P 47 
P 48 
P 49 
P 50 
P 51 
P 52 
P 53 
P 54 
P 55 
P 56 
P 57 
P 58 
P 59 
P 60 
P 61 
P 62 
P 63 
P 64 
P 65 
P 66 
P 67 
P 68 
P 69 
P 70 
P 71 
P 72 
P 73 
P 74 
P 75 
P 75 
P 76- 


SUtlROUTINE ERROR 

I ()1«0ATA( I I 
I l)2»OATA(2) 

ssn«o. 


WlUrt: (Otl5) 
I 


Q 1 
Q 2 
0 3 
Q 4 
0 5 
Q 6 
Q 7 
Q 8 
0 9 
0 10 
Q 11 
0 12 
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1 
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5 
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967 

6 
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C 
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980 

C 

981 

12 

982 

C 

983 

C 

984 


985 

C 

986 


987 


988 

c 

989 


990 

13 

99 1 
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C 

992 

15 

993 

16 

994 

17 

995 

18 

996 

19 

997 



IF TIME INCREMENTS NOT EQUALi INTERPOLATE 

IF (DT(IOH-OTI 1021) 1,8,2 

L=ID1 

K=ID2 

GO TO 3 

L=ID2 

K=I01 

M=I£NO(L) 

TIO=DT(K) 

T10H=0. 

DO 6 1*2, H 
TIDH=TIDH*OTU) 

IF (TID-TIOH) 5,4,6 
J= J+1 

C.FSI J)=0CFSI1 
TT0=T1D+0TIK) \ 

GO TO 6 
J = J+1 

^CFS( J ) *OCFS n-1 ,L) V( ( TID-TIOH+DT (L) ) / DT(L) ) *(OCFS ( I , L )-nCFS I I-l ,L) 

TIC=TID+DT{K) 

CONTINUE 

IEND(L)=J 

OT(L)*DT(K) 

DO 7 1*2, J 

OCFSI i,L) = CFS(n 

IF { lENOI lOll-IENDl 1D2) ) 9,9,10 

M* IF.NOI IDl ) 

GO TO 11 
M=1ENO(I02) 

T2*TIME 

DETERMINE ERROR 
DO 12 1=1, M 

ERR=OCFS( I ,ID1)-0CFS( I , ID2) 

WRITE (6,16) T2 ,OCFS( I , IDl ) ,OCFS( I , 102) , ERR 
T2»T2+0T{ IDl) 

SUM OF SQUARES OF ERROR 

SSE*SS6+ERR*ERR 

XM=M 

ERROR VARIANCE 
EVAR=SSE/XM 

ERROR STANDARD DEVIATION 
ESOEV’-SORTIEVAR) 

WRITE (6,17) ESOEV 

PERCENT ERROR FOR PEAK FLOWS 

ERPK=ABS( PEAK (IDl )-PEAK( 1021 ) 

PCTER=(ERPK/PEAK(IOl) 1*100. 

WRITE (6,18) PCTER 
PUNCH CODE 
IF (NPU) 14,14,13 
WRITE (7,19) ID1,ID2 
RETURN 

FORMAT! IHO ,T33, 'TIME* ,T55 j'FLOW 1*,T76,*FL0W 2*,T95,'ERR0R'/T34, 
1*HRS',T57,*CFS',T78,'CFS',T97,'CFS*) 

FORMAT (20X,F20.3t3F20*0) 

FORMATt IHO, no, 'ERROR STANDARD DEVIATION ^ SF10*3) 

F0RMAT(T10,*PEAK DISCHARGE ERROR == PERCENT*///) 

FORMAT! 'ERROR ANAL YS I S * , T2 1 , ' I D l=SIlfT29,MD 11=*, 111 
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998 


SUBROUTINE SECT 

R 

1 


C 

THIS SUBROUTINE COMPUTES THE SEDIMENT YIELD FOR A FLOOD 

R 

2 

999 


COMMQN CFSOOOI ,0CFS(300,6I , I END! 6 ) ,OATA ( 3 10) 1 0A( 6 » t DP ( 20 ) , NPU ,NHO 

R 

3 



1, SCFS120),CI20) ,A(20,6) ,Q 1 20 1 6 ) , RAI N ( 200 ) , DEEP ( 20 , 6) tNER i MAXNO.NCO 

R 

4 



2MM,ICC,NC0DE,DIST(6),SEGN(6) tCTBLEISOfll) ,IT6LEI50,2) |ZALFAI20),DT 

R 

5 



3( 6) f TI MEfPEAK(6 ) f ROINi ISG(6) 

R 

6 

1000 


I 0=DATA( 1 ) 

R 

7 

1001 


SOIL = DATAt 2) 

R 

8 

1002 


CROP=DATA( 3) 

R 

9 

1003 


CP=DATA(4) 

R 

10 

1004 


SL=OATA( 5) 

R 

11 


C 

COMPUTE SEDIMENT YIELD 

R 

12 

1005 


X=ROIN*DA( to) «53.333*PEAK(I0) 

R 

13 

100 6 


SED=95 56*S0I L*CRaP*CP*SL 

R 

14 

1007 


WRITE (6t3) SED 

R 

15 


c 

PUNCH CODE 

R 

16 

1008 


IF (NPU) 2«2il 

R 

17 

1009 

1 

WRITE (7,4) ID, SOIL. CROP, CP, SL 

R 

18 

1010 

2 

R E TURN 

R 

19 


C 


R 

20 

1011 

3 

FORMAT (lOX, 'SEDIMENT YIELD = FlO.l, • TONS') 

R 

21 

1012 

4 

FORMAT! 'SEDIMENT Y I ELD ' , T21 , • ID»' , 11 ,T29, • SOI L= ' , F5.3, T42 , • CROP 

R 

22 



!=• ,F5.3,T57,'CP = ' ,F5.3,T70, 'LS*' ,F5.3) 

R 

23 

1013 


END 

R 

24- 


//tOATA 
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ZALFA = 123A5678g0 - 


COMMAND TABLE 


START 12 

STORE HYD 2310 

RECALL HYO 3310 

COMPUTE HYO 4310 

PRINT HYO 5 2 

PUNCH HYD 6 1 

PLOT HYD 7 2 

AOD HYD 8 4 

STORE RATING CURVE 9100 

COMPUTE RATING CURVE 10310 

STORE TRAVEL TIME 11100 

COMPUTE TRAVEL TIME 12 5 

ROUTE 13 4 

ROUTE RESERVOIR 14100 

ERROR ANALYSIS 15 2 

sediment yield 16 5 

FINISH 17 0 
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